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We study neutrino-nucleon scattering and absorption in a dense, magnetized nuclear medium. 
These are the most important sources of neutrino opacity governing the cooling of a proto-neutron 
star in the first tens of seconds after its formation. Because the weak interaction is parity violating, 
the absorption and scattering cross-sections depend asymmetrically on the directions of the neutrino 
momenta with respect to the magnetic field. We develop the moment formalism of neutrino transport 
OQ ', m the presence of such asymmetric opacities and derive explicit expressions for the neutrino flux and 

^\ ■ other angular moments of the Boltzmann transport equation. For a given neutrino species, there is a 

ON ' drift flux of neutrinos along the magnetic field in addition to the usual diffusive flux. This drift flux 

depends on the deviation of the neutrino distribution function from thermal equilibrium. Hence, 
despite the fact that the neutrino cross-sections are asymmetric throughout the star, asymmetric 
Q ' neutrino flux can be generated only in the outer region of the proto-neutron star where the neutrino 

distribution deviates significantly from thermal equilibrium. The deviation from equilibrium is 
similarly altered by the asymmetric scattering and absorption, although its magnitude will still 
be quite small in the interior of the star. We clarify two reasons why previous studies have led 
to misleading results. First, inelasticity must be included in the phase space integrals in order to 
satisfy detail balance. Second, nucleon recoil must be included in order to find the leading order 
asymmetric cross sections correctly, even though it can be ignored to leading order to get the zero 
field opacities. In addition to the asymmetric absorption opacity arising from nucleon polarization, 
we find the contribution of the electron (or positron) ground state Landau level. For neutrinos of 
^/-^ , energy less than a few times the temperature, this is the dominant source of asymmetric opacity. 

Lastly, we discuss the implication of our result to the origin of pulsar kicks: in order to generate 
kick velocity of a few hundred km s _1 from asymmetric neutrino emission using the parity violation 
effect, the proto-neutron star must have a dipole magnetic field of at least 10 15 — 10 16 G. 
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■ I. INTRODUCTION 
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A. Astrophysical Motivation 

X, 

Neutrinos play an essential role in core-collapse supernovae and the formation of neutron stars GL^] . It is through 
neutrino emission that a hot proto-neutron star releases its gravitational binding energy and cools. The explosion itself 
also relies on the neutrinos for its success. The most important ingredient of neutrino transport in proto-neutron 
stars is the neutrino opacity in a dense nuclear medium. Much effort has been devoted to understanding various 
effects of neutrino-matter interactions at supra- nuclear density (e.g., |p|-|lT| and references therein). Neutron stars, 
however, possess strong magnetic fields. While the present-day, dipolar magnetic fields of most radio pulsars lie in 
the range of I0 12 — 10 13 G, it has been suggested that fields of I0 15 G or larger can be generated by dynamo process 
in proto-neutron stars |Lq ]. Several recent observations [p^-p"5| have lent support to the idea that soft gamma-ray 
repeaters and slowly spinning X-ray pulsars ("anomalous x-ray pulsars") in supernova remnants are neutron stars 
endowed with superstrong magnetic fields B > I0 14 G |H^,[l7|. It is therefore necessary to understand how neutrino 
opacities are modified by the presence of a strong magnetic field. This is the subject of our paper. 

A direct motivation of our study is to explore whether strong magnetic fields can induce asymmetric neutrino 
emission from proto-neutron stars to explain pulsar "kicks" . It has long been recognized that neutron stars have 
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space velocities that are about an order of magnitude greater than their progenitors' (e.g., [p^|Jl£ l). Recent studies of 
pulsar proper motion give 200 — 500 km s _1 as the mean 3D velocity of neutron stars at birth |20 |23[|, with possibly a 
significant population having velocity of order or greater than 700 km s _1 . Direct evidence for pulsars with velocities 
<; 1000 km s _1 comes from observations of the bow shock produced by PSR B2224+65 in the interstellar medium 
p4j , and studies of pulsar-supernova remnant associations [p5{ . A natural explanation for such high velocities is 
that supernova explosions are asymmetric, and provide kicks to nascent neutron stars. Support for supernova kicks 
has come from the detections of geodetic precession in the binary pulsar PSR 1913+16 |^,^| and the orbital plane 
precession in the PSR J0045-7319/B star binary and its fast orbital deca y |28| J29|| . In addition, evolutionary studies 
of the neutron star binary population imply the existence of pulsar kicks J3 C| — 133|] . 

Two classes of mechanisms for the natal kicks have been suggested. The first class relies on hydrodynamical 



instabilities in the collapsed stellar core j34 38 that lead to asymmetric matter ejection and/or asymmetric neutrino 
emission; but numerical simulations indicate that these instabilities are not adequate to account for kick velocities 
^> 100 km s _1 p^ , |39j,^0l . G lobal asymmetric perturbations in the presupernova cores are required in order to produce 
the observed kicks | [39|]41| . In this paper we are concerned with the second class of models in which the pulsar 
kicks arise from magnetic field induced asymmetry in neutrino emissions from proto-neutron stars. The fractional 
asymmetry a in the radiated neutrino energy required to generate a kick velocity Vicick is 



a = 0.028 



M 



I AM, 



3 x 10 53 



erg 



E„ 



Vk 



ick 



1000 km s~ 



(1.1) 



where M is the mass of the neutron star and -Etot is the total neutrino energy radiated from the neutron star. Since 
99% of the neutron star binding energy (a few times 10 53 erg) is released in neutrinos, tapping the neutrino energy 
would appear to be an efficient means to kick the newly-formed neutron star. Magnetic fields are naturally invoked 
to break the spherical symmetry in neutrino emission, but the actual mechanism is unclear. We first review previous 
work related to this subject. 



B. Review of Previous Work 

Beta decay in a strong magnetic field was first investigated in Refs. Jl2|,[l3| (See also Refs. Q,^5|). A number 
of authors have noted that parity violation in weak interactions may lead to asymmetric neutrino emission from 
proto-neutron stars p^|49|] . Chugai Q and Vilenkin [Q (see also Ref. Q ) considered neutrino-electron scattering 
and concluded that the effect is extremely smalfj] (e.g., to obtain Vkick = 300 km s _1 would require a magnetic 
field of at least 10 16 G). However, neutrino-electron scattering is less important than neutrino- nucleon scattering in 
determining the characteristics of neutrino transport in proto-neutron stars (e.g., |]5l| , p2| ). Similarly, Dorofeev et al. 
p7| considered neutrino emission by Urea processes in strong magnetic fields. However, as we shall see below (see 
Rcf . [53] ] ) , the asymmetry in neutrino emission is cancelled by the asymmetry associated with neutrino absorption for 
young proto-neutron stars where the neutrinos are nearly in thermal equilibrium. The size of the asymmetric flux 
due to absorption/emission processes is then dependent on the deviations from thermal equilibrium at the neutrino 
photosphere. 

Horowitz & Li j54j] suggested that large asymmetries in the neutrino flux could result from the cumulative effect of 
multiple scatterings of neutrinos by nucleons which are slightly polarized by the magnetic field. In particular, they 
found that the size of the asymmetry was proportional to the optical depth of the star to neutrinos (r ~ 10 4 ). The 
result was that kick velocities of a few hundred km s -1 could be generated by field strengths of only a few times 
10 12 G. Initial neutrino cooling calculation j5^| of a proto-neutron star in magnetic fields appeared to indicate that 
a dipole field of order 10 14 G is needed to produce a kick velocity of 200 km s _1 . The larger magnetic field required 
results from cancellations of the asymmetries associated with , v T and their antiparticles as well as the opposite 
signs of polarizations of neutrons and protons. Preliminary numerical study reported in Ref. [ [40"| drew a similar 
conclusion although it was claimed that only 10 13 G is needed to produce 200 km s" 1 . 

As appealing as the cumulative effect may be, we now believe that it does not work in the bulk interior of the 
neutron star [p5| , p6| , and the conclusions reached in Refs. ^,[53|,Q are incorrect J57|. In spite of the fact that the 
scattering cross-section is asymmetric with respect to the magnetic field for individual neutrinos, detailed balance 
requires that no asymmetric neutrino flux can arise in the stellar interior where neutrinos are in thermal equilibrium 



1 Note that Chugai's estimate for the electron polarization in the relativistic and degenerate regime (the relevant physical 
regime) is incorrect. This error leads to an overestimate of the effect as compared to Vilenkin's result. 
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to a good approximation. Since this issue is somewhat subtle and counter-intuitive, we discuss it in detail in §11 where 
we also point out where previous studies went wrong. 

A related, but different kick mechanism relies on the asymmetric magnetic field distribution in proto-neutron stars 
|^8|-]60|. Since the cross section for v e (D e ) absorption on neutrons (protons) depends on the local magnetic field 
strength due to the quantization of energy levels for the e" (e + ) produced in the final state, the local neutrino fluxes 
emerged from different regions of the stellar surface are different. It was found that to generate a kick velocity of 
~ 300 km s _1 using this mechanism alone would require that the difference in the field strengths at the two opposite 
poles of the star be at least 10 16 G. 

There have also been several interesting ideas on pulsar kicks which rely on nonstandard neutrino physics. It was 
suggested pH that asymmetric v T emission could result from the Mikheyev-Smirnov-Wolfenstein flavor transformation 
between v T and v e inside a magnetized proto-neutron star because a magnetic field changes the resonance condition 
for the flavor transformation. Another similar idea 62 relies on both the neutrino mass and the neutrino magnetic 
moment to facilitate the flavor transformation. More detailed analysis |6^,Q, however, indicates that even with 
favorable neutrino parameters (such as mass and magnetic moment) for neutrino oscillation, the induced pulsar kick 
is much smaller than previously estimated. We will not consider the issues related to nonstandard physics in this 
paper. 

Finally, we mention that previous calculations of neutrino processes in magnetic fields have generally neglected 
nucleon recoils (e.g., p0|j4^^ , ^7| , p3| , p4t ) . Although this simplification is justified in most cases without magnetic 
field, it is invalid in a magnetic field because the asymmetric part of the opacity depends sensitively on the phase 
space available in the scattering/absorption. This and other technical issues (such as Landau levels) will be addressed 
in our paper. 



C. Plan of This Paper 



In this paper, we carry out a systematic study of neutrino-nucleon (i>, N) scattering and electron neutrino absorp- 
tion/emission (y e + n ^ p + e~ and D e + p ^ n + e + ) in a dense, magnetized nuclear medium. These are the most 
important sources of opacity for neutrino cooling of the proto-neutron star in the first tens of seconds after its forma- 
tion, when most of the binding energy of the neutron star is radiated as neutrinos. Our study goes beyond merely 
calculating differential cross-sections of the neutrinos in magnetized medium in that we derive the expression for the 
neutrino flux and other angular moments from Boltzmann equation. This is necessary in order to determine whether 
the effect of parity violation introduces any asymmetric "drift flux" in addition to the usual "diffusive flux" . Indeed, 
there are a number of subtleties in these derivations that several previous studies have arrived at incorrect results (see 
§I.B and §11). We show that, despite of the fact that the scattering/absorption cross-sections are asymmetric with 
respect to the magnetic field for individual neutrinos, there is no "drift flux" when the neutrinos are in thermal equi- 
librium; the drift flux is proportional to the deviations from equilibrium, which are small below the neutrino-matter 
decoupling layer. Hence, asymmetric neutrino emission can be generated only near the surface layer of the star. 

In §11 we discuss a simplified calculation and point out its problems. This serves as an illustration of various issues 
that one must pay attention to in order to obtain the correct answers. We begin our formal theoretical development 
in §111, where the relevant cross sections are defined starting from the Boltzmann transport equation. It is important 
to distinguish different cross sections (e.g., those related to scattering into the beam and scattering out of the beam in 
the Boltzmann equation) in order to satisfy the principle of detailed balance, which states that in complete thermal 
equilibrium the collisional term in the Boltzmann equation vanishes. 

Section IV contains a detailed calculation of v — N scattering in magnetic fields. Starting from the weak interaction 
Hamiltonian, we compute the scattering opacity, carefully including the effect of nuclear motion to lowest nonvanishing 
order. This opacity is then used to find the contribution to the angular moments of the transport equation. Explicit 
expressions are obtained for the outer layer of the neutron star where nucleons are nondegenerate and where parity- 
violating asymmetric flux can be generated. A technical complication arises from the effect of small inelasticity: Even 
in the regime where the nucleon recoil energy is much smaller than the neutrino energy and temperature, phase space 
considerations require that the inelasticity effect be included when deriving the asymmetric flux and other moments 
(The situation is similar to that found in the derivation of the Kompaneets equation in electron-photon scattering; 
see Ref. ||). 

In Section V we calculate the cross sections for the absorption of v e and v e by nucleons. As in the scattering case, it 
is necessary to include nucleon recoil in the absorption calculation. Additional complications arise from the quantized 
Landau levels of the final state electrons (or positrons). In particular, the ground state Landau level of the electron 
introduces an effective electron "polarization" term in the asymmetric part of the opacity. For certain parameter 
regimes, this electron "polarization" term dominates the asymmetry in the absorption opacity. We also demonstrate 
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explicitly that the Landau levels of protons have no effect on the absorption opacity since many levels are summed 
over for the situation of interest. 

In Section VI we combine the results of §IV and §V to derive the angular moments of the Boltzmann transport 
equation. These moment equations are truncated at the quadrupole order, since we expect that the contributions of 
the higher order terms to the asymmetric flux are smaller. As expected, our explicit expression for the neutrino flux 
contains the usual diffusive flux plus a drift flux along the magnetic field. This drift flux, however, depends on the 
deviation for the neutrino distribution from thermal equilibrium. Although strictly speaking our truncated moment 
equations break down near neutron star surface, these equations are accurate below the neutrinosphere, and provide a 
reasonable physical description of the neutrino radiation field throughout the star. Finally in §VII we use the moment 
equations to obtain an order-of-magnitude estimate of the asymmetric neutrino emission from the proto-neutron star 
due to the parity violating processes considered in this paper. 

Throughout this paper, we treat nucleons as noninteracting particles This is clearly a simplifying assumption. 
In reality, strong interaction correlations may significantly change the neutrino opacities (e.g., |lC| , |lT[ ] and references 
therein). However, the goal of this paper is to consider whether there is any new effect associated with strong magnetic 
fields. For this purpose, it is certainly appropriate to focus on noninteracting nucleons, particularly since there are still 
large uncertainties in our understanding of matter at super-nuclear densities. Moreover, for certain nuclear potentials, 
the nuclear medium effects merely amount to giving the nucleon an effective mass, and therefore the result of this 
paper can be easily extended. For general nuclear interactions, it is likely that the qualitative conclusion reached in 
this paper will remain valid, although this issue is beyond the scope of this paper. 

Unless noted otherwise, we shall use units in which h, c and the Boltzmann constant ks are unity. 



II. A SIMPLIFIED CALCULATION OF THE SCATTERING EFFECT AND ITS PROBLEMS 

As indicated above, the effect of asymmetric neutrino-nucleon scattering is sufficiently subtle that several previous 
studies have led to a wrong conclusion. It is therefore instructive to consider a simplified treatment of the problem to 
understand how previous work went wrong. This section also serves as an illustration of the various issues that one 
must pay attention to in doing such a calculation (Our systematic calculation is presented in §§III-VI). 



A. The Calculation 



We shall follow the treatment as given in Ref. [|53[ (Refs. [|40 54 used a Monte-Carlo method for the neutrino 
transport, which is less transparent for our analysis). For the purpose of clarity, we consider the scattering of 
neutrinos by nondegenerate neutrons. Assuming the scattering is elastic (which is a good approximation since the 
neutrino energy k is much less than the neutron mass), one can easily obtain the matrix elements when the initial 
neutron has spin along the magnetic axis (the z-axis): 



cosi(0'-0) + Acosi( 



l^kT^k'tl — 2G F Cy 

l^kT^J 2 =2G 2 F c\ (2 sin cos ±e) 



(2.1) 
(2.2) 



where k (k') is the initial (final) neutrino momentum, 9 (9') is the angle between k (k') and the z-axis (assuming 
azimuthal angle 4> — 0), and Gf , cy, ca, A are weak interaction constants as defined in Appendix A. The differential 
cross section, for the initial nucleon with spin along the z-axis, is given by 



da 



k 2 



(2tt) 2 
fc 2 

(2tt) 2 



l^kT^k'tl 2 + l^kt^k'il 2 



r 2 r 2 



J F L -V 



(1 + 3A 2 ) + 2A(A + 1) cos 9 - 2A(A - 1) cos 9' 
(1- A 2 )cos(6>-6»') 



(2.3) 



One can similarly obtain (c?cr/dr2')kx^k' when the initial nucleon spin is anti-parallel to the z-axis. For general nucleon 
spin polarization P — (cr z ), the differential cross section is given by 
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da 
dtt' 



1 + P\ f da 

~dn* 



kT^k' 



1 - P\ f da 
dM 



kl^k' 



(1 + 3A 2 ) 



1 - A 2 

1 + e in cos 9 + e out cos 6' + — -r cos(# - 9') 

1 + 3\ z 



where 



= 2P 



A(A + 1) 
(1 + 3A 2 ) : 



Eout — — 2P 



A(A-l) 
(1 + 3A 2 )' 



(2.4) 



(2.5) 



This clearly indicates that the scattering is asymmetric with respect to the magnetic field, a direct consequence of 
parity violation in weak interaction. A similar expression was derived in Ref . |54| , although a different sign in ei„ and 
e ut was given. 

Next we study the consequence of the asymmetric cross section on the neutrino flux. The Boltzman transport 
equation for the neutrino distribution function /^(k) can be written in the form: 



dt 



n-v/„(k) = /dn'(^) t f*W[i-f*Q4] 



dn 
dti/ 



/„(k) [l-/„(k')] 



(2.6) 



where SI and SI' are unit vectors along k and k', respectively, and the (elastic) differential cross section per unit 
volume can be written in the form 



dn 
dfl' 



n->n' 



K 
47T 



1 + einft ■ B + 

^out 

SI' B 



const. (SI ■ SI') 



(2.7) 



with B the unit vector along the magnetic field. Note that in eq. (2.6) we have neglected ne utrin o absorption for 
simplicity. The first order moment of the transport equation is obtained by multiplying eq. (2.6) by SI and then 
integrating over dil. The specific neutrino flux is then given by^ 



F„ = ~VU V + - (eout 

6K 6 



e iB )cU v (l- f u )B, 



(2.8) 



where U v is the specific neutrino energy density. According to eq. (2.8), the neutrino flux consists of the usual 
diffusive flux, P<jiff oc V£/„, and a "drift" flux P drift along the magnetic field. The drift flux induces asymmetric 
neutrino transport. One can easily see that the ratio -Fdiff /-Fdrift is of order (e ou t — £in)'?~, where r ~ kR is the optical 
depth of the star (R is the stellar radius), i.e., the asymmetry increases with r. This is the origin of the cumulative 
effect discussed in Refs. EolMRl. 



B. The Problems 

The calculation presented above, while physically motivated, is actually incorrect. There are two problems: first, 
the asymmetry terms (those proportional to e; n and e ou t) in eq. (|]J) are incomplete. Even in the regime where the 
elastic approximation is highly accurate from the energy point of view, small inelasticity will affect the asymmetric 
part of the cross section. This comes about because the asymmetric terms depend in a subtle way on the available 
phase space of the scattering. Indeed, our full calculation presented in §IV reveals additional terms in the expressions 
of £j n and e ut- Moreover, to obtain the the correct expression for the asymmetric neutrino flux, the elastic cross 
section is inadequate; one must incorporate the full inelastic effect in the Boltzmann equation (see §IV). A similar 
comment can be made for neutrino-nucleon absorption, where one must incorporate the recoil motion of the nucleon 
as well as the Landau levels for the electron in order to obtain the correct cross section (See §V). 



An overall factor of (1 — const. /3) 1 has been dropped in this equation. 
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Second, and more importantly, eq. (2.6) i s in correct, therefore eq. (2.8) is also incorrect and there is no drift flux 
proportional to U u , The problem with eq. (|2.6| ) can easily be seen by considering detailed balance (e.g., Ref. @): 
The right-hand-side of eq. (2.6) must vanish when neutrinos are in thermal equilibrium with the matter. Substituting 

/i/(k) and / y (k') by the equilibrium distribution f„°\k) and using eq. (2.7), we find the RHS of eq. (2.6) to be 

ftfv Jl — ji^Xeout — £in)J~2 • B. It is exactly this violation of detailed balance that gives rise to the drift flux term in 
eq. ([2^) . It is also clear that any physical drift flux must depend on the deviation from the equilibrium distribution. 

Equation ( |2.6| ) is the starting point of almost all astrophysical radiative transport studies (e.g., Ref. |^|). However, 
it is invalid in the presence of asymmetric scattering. A crucial (but incorrect) assumption implicit in eq. (|]^) is that 
the cross section for scattering into the beam (propagating along (dn/dtyw^si, is related to that for scattering 
into the beam, (dft/<if2')n— >n', by merely switching Jl and fi'. In reality, how ever , the two cross sections have slightly 
different forms such that detailed balance is satisfied in equilibrium (see eq. Q3.9[ ). In other words, although the first 
(second) term in eq. ( |2.6| ) represents a good approximation to the actual probability of scattering into (out of) the 
beam, the error in their difference is significant. It will become clear in our study presented the following sections 
that to properly take into account of the detailed balance constraint, one must incorporate inelasticity — no matter 
how small — -into the Boltzmann equation when deriving the asymmetric neutrino flux. 

III. GENERAL FORMALISM 

In this section, we set up the general framework to study neutrino transport in magnetic fields. As §11 shows, it is 
important to properly define the relevant cross sections which enter the transport equation. The actual calculations 
of the cross sections are given in §IV and §V. 

The Boltzmann equation for neutrino transport is written in the form 



dt 



n v/„(k) = 



\dMky 








dt 


+ 

sc 


dt 


abs 



(3.1) 



where k = kfl is the neutrino momentum, both scattering and absorption collisions terms are included on the right- 
hand-side of the equation and we have suppressed the position and time dependence in the neutrino distribution 
function /„. 



A. The Scattering Term 

The collision term due to neutrino- nucleon scattering can be written as (e.g., Ref. Ipql ) 



dt 



E 



d 3 k' d 3 p d 3 p' 



(2^) 4 <5 4 {P + K-P' - K 1 ) \M SS , (n, «')f 



(2tt) 3 (2tt) 3 (2tt) 3 



x [(1 - /„)(! - Mf'ttfl - fuf N {l - f' N )(l - f v ) 



(3.2) 



where f v = f v Qs) and f v = /j/(k') are the initial and final state neutrino distribution functions, /jy = In(E) and 
/at = Jn(E') are the initial and final state nucleon distribution functions, P (P 1 ) and K (K 1 ) are the initial (final) 
state nucleon and neutrino four vectors, respectively, and s, s' — ±1 specify the initial and final nucleon spins. Time- 
reversal symmetry has been used to relate the matrix element for scattering in and out of the beam. Note that in the 
case of neutrino-proton scattering, we neglect the Landau levels of proton, and therefore the proton momentum is a 
well-defined quantity. This is justified because many Landau levels are occupied for the conditions in a proto-neutron 
star, and the change in the available phase space due to the Landau levels is negligible. Nucleons are always in thermal 
equilibrium, and the nucleon distribution function is given by 



ex P [{E- m )/T] + l 



(3.3) 



where /ijv is the nucleon chemical potential (excluding rest mass). As the neutrinos exchange energy with matter 
only through the weak interactions, their distribution can deviate from equilibrium, especially in the outer layer of 
the star. 

The scattering rate can be rewritten in a more conventionial form as follows. Define the differential cross section 
(per unit volume) to be 



6 



dT 
dk'dtt 1 



^Y,\M ss ,(n,n'f S ss ,(q ,q), 



where the "nucleon response function", S ss i (qo, q), is 



S aa > (go,?) 



d 3 p d 3 p' 



(2tt) 4 5 4 (P + K-P'- K') Ml - /^). 



(2tt) 3 (2tt) 3 

Here we have defined the energy transfer qo and the momentum transfer q via: 

q = k-k', q = |k - k'| = (fc 2 + fc' 2 - 2fcfc'fi • 0') 1/2 . 
For nucleons in thermal equilibrium, energy conservation gives 

(i - Mf' N 



fN(l-f' N ) 



= exp[(£ - £')/T] = cxp[-(fc - k')/T] = exp(-g /T). 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



Using this expression to relate the scattering into and out of the beam and plugging in the differential cross section, 
the scattering rate becomes 



Ot 



dk' IdVt' 



dT 



dk'dW 



(3.8) 



Note that one can also explicitly define a differential cross-section, [dF/(dk'd£l')]+, for scattering into the beam by 
writing 



df(k) 



dt 



dk' dfl' 



{dk'dn)^ 1 u)f '» {dk>dn>) M1 ! ' v) 



Clearly we have 



dT 
dk'dVL 



dT 



dk'dQ' 



(3.9) 



(3.10) 



All the microphysics is n ow c ontained in dT /(dk'dVt'). It depends only on k, fc', and S~2 ■ S~2' , or equivalently, fc, qo and 
q, as can be seen from eq. ( |3.5[ ). Note that if one sets qo/T — 0, then the neutrino Fermi blocking terms proportional 
to / 1/ (k)/ zy (k') cancel (e.g., Ref. |(||). Hence, the neutrino degeneracy does not enter the scattering rate if the elastic 
limit is taken in the phase space i nteg rals. 

It is instructive to compare eq. (3J;) with the (wrong) eq. ( |2.6| ). In general, the cross section for scattering into the 
beam has a different form as that for scattering out of the beam. This difference, even though numerically small, is 
essential for maintaining detailed balance in thermal equi libriu m (see §111. C below). Also, as we show in the next few 
sections, one cannot trivially take the elastic limit in eq. ( [3.10 ), because this will lead to zero drift flux even when the 
neutrino distribution deviates from thermal equilibrium. 



B. The Absorption Term 



The collision term in eq. (3.1) for absorption/emission is 



\dMky 




dt 


-J 

abs J 



[fpfe(l 



/n)(l-/,)-(l-/p)(l-/e)/,/n] 



(3.11) 



where W^ hs ^ is the transition rate (S-matrix squared divided by time) for absorption, and we have used time-reversal 
invariance of the weak Hamiltonian to relate the S-matrix for each direction. The notation dH denotes sum of states 
(including spins). Note that since we will include Landau levels for electrons and protons, dX\ e _ p is not equal to 
<^ 3 Pe,p/(27r) 3 (see Appendix C). The components of the transverse momentum (perpendicular to the magnetic field) 
are not conserved, although we still have conservation of the z-momentum and energy conservation 



k + E n + Q=E e + E. 



(3.12) 
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where Q is the mass difference between neutron and proton (Recall that E n and E p do not include rest-mass). Since 
the electron, proton, and neutron are in thermal equilibrium with Fermi-Dirac distributions, we have the equality 



/p/e(l fn) 



exp 



k — /i„ 



(l-/p)(l-/e)/„ 

where we have defined the neutrino chemical potential 

IX V = H e + [l v - fJL n - Q. 



T 



Equation (3.11) then takes on the standard form 



Of 



K< ahs ^6f u , 



(3.13) 



(3.14) 



(3.15) 



a.bs 



where 8f v = /r/(k) — fv (k) measures the deviation of neutrino distribution from thermal equilibrium (see below). 
Here 



K *(abs) _ ^(abs) 



1 + exp 



Hv - k 
T 



and 

K (abs) ig 

the absorption opacity: 



,(abs) = / d n e dn nC ?n p iy/; bs) (i - / p )(i - / e )/„ 



The factor [1 + e^*" fe )/ T ] in K*( abs ) takes into account the effect of stimulated absorption (e.g., Ref. |7Q]). 



C. Detailed Balance 



(3.16) 



(3.17) 



In thermal equilibrium, the neutrino has the Fermi-Dirac distribution function 



fu = fi°\k) = 



1 



We then find 



/5 0) (i-CT 
(i-/i 0) )/; (0) 



/(Oh 



expp-^)/T] + 1 



exp[(fc' - A)/T] = exp(-<to/r), 



so that the scat tering rate in cq. ( [3.8] ) is zero, as required by detailed balance. Similarly, for /„ = fj?\ the absorption 
rate in eq. (3.15) vanishes. Therefore the only nonzero contribution to {df v /dt) sc and (df^/dt)^ must be proportional 
to the deviation of neutrino distribution from thermal equilibrium. This implies that there is no drift flux along the 

magnetic field proportional to (see §11. A) . 



(3.18) 



(3.19) 



As noted before (§111. A), when one writes the scattering rate in the form of eq. (3.9), it is essential to distinguish 
(dT I dk' 'dft) + from (dT/dk'dDf) in order to satisfy detailed balance. 



D. Deviation from Thermal Equilibrium 



In order to calculate the size of the asymmetric flux, we must consider the deviation of neutrino distribution from 
thermal equilibrium: 



8f u (k) = /„(k) - 



(3.20) 



For v e and i> ei the neutrino-matter energy exchange is mediated primarily by absorption and emission via the URCA 
processes (and to a lesser extent by electron- neutrino scattering), while neutrino transport is affected by both absorp- 
tion and v — N scattering. The result is that the decoupling sphere of electron type neutrinos lies only slightly deeper 
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than the neutrinosphere and only a small region over which the cumulative asymmetry can develop. For /i and r 
neutrinos, the transport opacity is primarily from v — N scattering, while energy exchange is due to inelastic v — e~ 
scattering. As a consequence, the decoupling layer is much deeper than the neutrinosphere and the asymmetric flux 
has a large optical depth over which to develop. Unfortunately, as we will show, the flux asymmetry due to the /i and 
r neutrinos is cancelled by that from the corresponding antincutrinos. 

For the purpose of deriving the moment equations of neutrino transport, we shall expand 5/„(k) in spherical 
harmonics up to quadrupolc order as 



<y/„(k) = g (k) + 3n ■ h(k) + y/y(fc)P y (n) 



(3.21) 



where 



Vi 



(3.22) 



(Th e com ponents of Vij can be explicitly expressed solely in terms of the quadrupole spherical harmonics Y2 m ). In 
cq. ( 3.21 ), g(k) is the spherically symmetric deviation from equilibrium, h(fc) represents the dipole deviation which 
leads to the flux, and Iij(k) is the tensor describing the pressure asymmetry. Since Vij is symmetric we can choose 
Iij to be symmetric, leaving six independent elements. Moreover, we shall choose iy to be traceless (The nonzero 
trace can always be incorporated into g v ). 

To be explicit about the physical meaning of each component of £/„, one can relate g, h, and iy to more commonly 
used quantities. The energy density per unit energy interval is 



k 2 dn 



rkfu 



4irk 3 



ti 0) (k)+g(k) 



(3.23) 



The energy flux per unit energy interval is 



, k 2 dn,„ p 471-fc 3 , , IN 
F fe = / TT^k^lU = TT^Hk)- 



(2tt) 3 



(2tt) 3 



(3.24) 



Using the indentity 



—UilljPki = -—Otjdki + — (dikdij + diidjk) 



(3.25) 



and the fact that 7y is traceless, the pressure tensor per unit energy interval is 



k 2 dn , n , _ 47rfc 3 

kViQjU - j—^ 

2 47Tfc 3 .,. 



(27T) 3 



3% 



(4 0) (fc)+ 5 (fc) 



(3.26) 



Thus Iij is the anisotropic portion of the pressure tensor. 

Note that we have truncated our expansion at the quadrupole order angular dependence, since each successive term 
will be smaller than the previous by a factor of |nj ~ r _1 (where r is the optical depth). It will be shown that 
_1 , and g and Z$j both scale as fl t~ 2 (for B = 0). The ?-th spherical harmonic would have coefficients 



h oc 



which scale as fi, r 



IV. NEUTRINO-NUCLEON SCATTERING 



In this section, we calculate the differential cross-section for v— N scattering in magnetic fields for general conditions 
of nucleons (non-relativistic but arbitrary degeneracy) . We also obtain explicit expressions for the angular moments 
of the scattering term of the Boltzmann equation. These moments are then used in (§VI) to obtain the neutrino flux, 
as well as the spherical and quadrupole deviations from thermodynamic equilibrium, on which the asymmetric flux 
depends. 

We approximate the nucleons as nonrelativistic particles with energy (excluding rest mass) 
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E(p,s) = j^- Sf iBB, (4.1) 

where B is the magnetic field strength, s = ±1 is the nucleon spin projection along the z-axis (in the direction B), 
and [Ib — geh/ (2mc) is the nucleon magnetic moment {g n = —1.913 for the neutron and g p = 2.793 for proton). In 
this section, we shall omit the label n or p whenever possible, denoting final-state quantities by a prime. Also, the 
quantization of the proton energy levels and nucleon- nucleon interactions are neglected (see §11. A) , although we shall 
include proton Landau levels in our calculation of the absorption opacity (§V). 



A. The Differential Cross Section 



The differential cross section, defined in eq. (3.4), can be evaluated analytically to linear order in B for general 
conditions of nucleons. Since /i^B = 3.15 x 10~ 4 gi?i4 MeV (where B^ is the field strength in units of 10 14 G) is 
much smaller than the temperature or nucleon Fermi energy in the proto-neutron star, an expansion in the lowest 
nonvanishing power of B is an excellent approximation. 

The matrix element, \M sa '(ft, J~i')| 2 , for the ca se in which both the initial and final nucleon states are polarized has 
been derived in appendix A with the result (eq. [ All ): 



\M ss ,(ft, ft')f = ^g 2 f c 2 v { (i + 3A 2 ) + (i - a 2 ) ft • ft' 



-2A(A + l){sft + s'ft') ■ B - 2A(A - l)(sfl' + s'ft) ■ B 
(1 - A 2 ) (1 + ft ■ ft') + 4\ 2 ft ■ BO' • B | 



(4.2) 



where Gf, cv,and A = ca /cv are the weak interaction constants defined in appendix A. The nucleon response 
function, defined in eq. ( |3.5| ), has been calculated in Appendix B. It can be written as S ss ' = Sq + 5S ss i, where So is 
the spin-independ ent B = result and SS SS ' is the correction linear in B. Combining the expressions for \M ss (ft, ft')\ 2 
and S ss i into eq. (3.4), we find: 



dT 



— = A (k, k', ft') + SA+ (k, k', n') ft ■ B + 6A- (fc, k', p') ft' ■ B, 



dk'dVL 



(4.3) 



where fi' = ft ■ ft' (no t to be confused with the nucleon magnetic moment, ps, or the nucleon chemical potential /xjv) 
The first term in eq. (4.3) is the B = result: 



<„/2 



A (MV) 



(2tt)= 



J2\M ss ,(fl,ft')\ 2 S (q ,q) 



:2G 2 F c 2 v [(1 + 3A 2 ) + (1 - A 2 ) M '] S (qo, q), 



(4.4) 



with 



So(qa,q) 



m 2 T 1 
2irq 1 — e~ 



In 



1 +exp(-x ) 



1 + exp(— xq — z) 



(4.5) 



and we have defined 



x 



(g - q 2 /2m) 2 _ fj, N 
AT(q 2 /2m) T 



and 



Qo 
T ' 



(4.6) 



The second and third terms in eq. (4.3) correspond to the corrections arising from nonzero B: 



6A + ft ■ B + SA^ft' ■ B = ^3 \Mss>(fl, fl')\ 2 5S ss ,(q , q), 



(4.7) 



with the coefficients: 
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5A± {k, k', n') 



k' z 2G 2 F c 2 v m 2 X^ B B 



(27T) 



jrq 



[exp(a;o) + 1] [1 + exp(-a; - z)] 



1 ± A 



<Z 2 y 



(4.8) 



The reason for writing the cross section in the form of eq. (4.3) is that the angul ar de pendence needed to find the 
moment equations (see §IV.C) is now manifest. Note that the cross section in eq. (4.3) exhibits parity violation. If 



the parity operation is taken, the vectors £1 and ft' reverse sign and the pseudovector B keeps the same sign so that 
the cross section does not retain the same form. Also note that the cross section for scattering from the state to 
the state SI' does not have the same numerical value as the reverse process. However, this does not mean that time 
reversal invariance is violated. The inequality arises from averaging the matrix element over the nucleoli distribution 
functions. Indeed, the matrix element in eq. (4.2) can be explicitly shown to satisfy time reversal invariance by 
simultaneously interchanging all initial and final state labels. 



B. The Differential Cross Section: Nondegenerate Nucleon Limit 



Even after expanding the cross section in Eqs. (4.4) and (|h§|) for small magnetic fields, the expressions are still 
quite difficult to evaluate in general. However, as discussed in §111, asymmetric drift flux can develop only when 
the neutrino distribution deviates from thermal equilibrium (i.e., above the decoupling sphere). This occurs in the 
regime where nucleons are nondegenerate (at density p ~ 10 12 — 10 13 g cm -3 ). In this subsection we derive simplified 
expressions of Aq and SA± which will be useful for obtaining the angular moments of the scattering term (§IV.C and 
D) and neutrino flux. 

For degenerate nucleons, the characteristic neutrino energy transfer in each scattering is of order qo ~ fc(T/m) 1 / 2 <C 
k. The cross section peaks sharply around k' = k, and we can evaluate Aq, SA± in a series in the small parameter 
(T/m) 1 / 2 . Define the dimensionless quantities 



e = [4(1 - n')T/m] 



1/2 



k 



ck 



(4.9) 



so that the range of u, the dimensionless neutrino energy, is from — 1/e to oo. Using the expansion of the nucleon 
response function derived in Appendix B, we have, to linear order in e, 



-1m ^ [^^J [1 + 3A 2 + (1-AV] 



ken 1 / 2 



1 



and 



Gpcyk\ 2[i B B n 

~f~ k^ 6 
1 , 



2tt 
Aem 

k{l-n')\ 2 



1 + —eu + eu 



3 2 



fc(l-M'). 



fc(l-//) 



em 



(4.10) 



,4 _ 2 M m2 

em 



In deriving these expressions we have used the the B = equation 

2^ 3 n 1/2 



exp 



i 3 T 3 



(4.11) 



(4.12) 



to relate the nucleon chemical potential /ijy to its number density n (The corrections due to finite B are of order B 2 ). 
These expansions of Aq, 5A± are valid under the conditions (see Appendix B for details) T -C m, k -C (mT) 1 / 2 , 
HbB <C T and k > k mm = /iei^m/T) 1 / 2 ~ lO~ 2 \g\Bi±T~ 1 / 2 MeV. These conditions are satisfied for the conditions 
of interest in our study^J. 

Note that dimensionally, 5A± is smaller than Aq by a factor of order (ibB/T, but also note that the quantity em/k 
(which appears on the second line of eq. [4.11]) is of order V mT /k ~ Wrn/T, which can be quite large. This point 



will be important when we consider the size of the neutrino drift flux. 



3 For k < k ml 
made (the k 



n , the expansions leading to eqs. (4.10) and (4.11) are no longer valid, so that different approximations must be 
limit). Since this is a relatively small range of neutrino energy, we ignore this complication here. 
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C. Moments of the Scattering Rate 



To derive the expression for the neutrino flux, one needs to take the angular moments of the Boltzmann transport 
equation (i.e., multiply the equation by some power of $7 and then integrate over d£l). In this subsection we derive 
the general expressions for the moments of the scattering rate in the Boltzmann equation. In the next subsection we 
shall evaluate these expressions explicitly for the regime when the nucleons are nondegenerate and the scattering is 
approximately elastic. 



We first write the scattering rate (eq. [3JS]) in terms of 8f v = fv — fv ■ 



at 



dk' /dn'^j^y [C{k, k')6fl + D(k, k')SU + E(k, k')Sf'JU] 



(or 



2T 



where the C, D, E coefficients, and their qo/T <C 1 expansions, are 

C(k, k') = e-WT (j _ ; (o)) + f (0) ^ x + (_| + (1 _ /(o)) 

D(k, k') = [e-to^fW' + 1 - fl^'} * 1 + §/<°> - J, (f^ ■ 

E(k,k') = l-e-^ T ^^-^ r2 . 

The nonlinear terms E{k, k^Sf^Sfl in eq. ( 4.13[ ) will be dropped since we consider the regime where the deviation of 
/„ from thermal equil ibriu m is relatively small (the regime where the moment formalism is valid) . 

Plugging eqs. (3.21) and (4.3) into eq. ( 4.13Q and then performing the azimuthal integrals using the indentities 



dJl_ 
dk 



(4.13) 

(4.14) 
(4.15) 
(4.16) 



— ili = fi il t 

Z7T 



''^njn^ = p 2 (p')n l n J + p 2 (//)) 



2vr 



2tt 



fc<5y] 



(4.17) 



where the three-index tensor equivalent to the Y^ m is defined by Vijk{n) — (5fljfljfife — fliSjk — — f2fc&y)/2, we 

find 



5/,(k) 



<9t 



cifc' / d/i' 



A + SA+n ■ B + 6A_n' ■ B 
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C ( g' + 30' • h' + —4^,- ) + £> ( <? + 3H • h + —h 3 V n 



10 



2tt 



dfc' / d/i' 



A + SA+fl ■ B 



(C 9 ' + Dg) + 3 (C/x'h' + Dh)-n 



— [CP^Wij + Dh^ 



+2tt 



d/z'<L4_ { jtx' (C 5 ' + £> ff ) ft • B + 3CP 2 (/i')(^ ■ B)(« • h') 



+C [1 - P 2 (//)] B ■ h' + 3D/i'(n • B)(n • h) 



~Y c ^ lI ' jk 



3 1 3 

gPijfe - -^n'QiSjk + yr// + QkSij) 



(4.18) 
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To calculate the moments of the scattering rate, the following identities are needed: 



dn 
An' 



1 



dn 

~An~ Vl - 



dn 

An 

dn 

An 



VijVki 



n^nAVki 



1 3 

-zSijSki + — {SikSij + SuSjk) 
10 20 

- — SijS k i + — (Sikdij + 5ii8 jk ) 
15 10 



T~ V »k = 
An 

dn 

An 



'Pijk^l — 0. 



(4.19) 



Also note that any integral of an odd number of fij's over the solid angle gives zero. The zeroth moment of eq. (4.1 h 
is 



dn 

An 



at 



2n / dk' / dfi' A (Cg' + Dg) 



6A + (C(i'ti + Dh) ■ B + 5A_(Dn'h + Cti) ■ B 



(4.20) 



the first moment is 



An 1 



at 



pOO pi 

2n I dk 1 I dn'A [Cfi'h'i + Dhi] 
Jo J -i 



2n 



dk' / dp' {Cg' + Dg) (8A+ + p'SA^) B, 



a pOO pi 

— J dk' J ^ d^' [SA + (CP 2 (m')/« + DI i:i ) + n'SA- (CI'j + DI^)} B v (4.21) 



and the second moment is 



/ m = 2n r dk ' s\ dfi>A ° [cP2 w>ij + m ^ 



2tt— / dk' / du! 
10 Jo J-i * 



C (n'5A+ + P 2 (n')6A_) [ h^Bj + h'jBt - -<%h' • B 



D (SA+ + n'&AJ) ( hiBj + hjBi - -^h ■ B 



(4.22) 



D. Moments of the Scattering Rate: Nondegenerate Nucleon Limit 



We now evaluate the moments in eqs. ( 1.20h -( 4.22) ex plicitly for nondegenerate nucleons (near the stellar surface) 
using the e xpress ions of Ap and SA± (eqs. 4.10[ and [4.11]) as derived in §IV.B for small inelasticity. After substituting 
eqs. ( 4.14 ) and ( 4.15 ), it will be necessary to evaluate moments of qo/T against Aq and 8A±. These moments are 
defined as 



MS = dk'A (k, k', M ') (I) " = ek (-^j " f 



duAou r ' 



and 



5MI = 



dk'SA ± (I)" 



ek 



-ek 



du5A±u n \ 



(4.23) 



(4.24) 
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where we have used the eq. ( |4.9D and in the du integrals we have extended the lower limit (— e 1 ) to — oo (since e<l) 
Only the following values will be needed f^J 



Mn 



Ml 



Mt = 



^ G F Cyk J 
GpCyk 
GpCyk 

2tt 



n[l + 3A 2 + (l-A 2 )//l 



n[l + 3A 2 + (l-A 2 yl 



?i[l + 3A 2 + (l-A 2 )//] 



-(1-/0 



T 



k 2 

2-^(1 -M') 
ml 



8M\ 



/ G F c v k \ 2 2n B B 



V 2tt / 



T 



A ± A 2 T 4A 2 



,T 



± 



GfcvA:\ 2^ B i? 2 
n — — — 2A . 



2tt 



1 



T 



(4.25) 



Let the notation 0(n) mean a term which contains a factor of (qo/T) n in the integrands of eqs. ( 4.20] )-( 4.22| ). In 
the zeroth moment, the O(0) term is zero and Mq is of the same size as Mq so we need to expand to 0(2) for the 
B = piece. For the B ^ piece, the 0(0) term is not zero, but SM± is of the same size as 8M± so we need to 
expand to 0(1). The zeroth moment is then given by 



dft 
4tt 



dt 



k 



G 



dk 



dk [ 



dk 



4^h I B. 



(4.26) 



where we have defined the zero-field scattering opacity (per unit volume) 

2 



4 8C) = y(^) (l + 5A 2 )n, 



(4.27) 



and the dimcnsionless asymmetry parameter 



6A 2 



Ms B 



(1 + 5A 2 ) T 



(4.28) 



For the first moment, the B = terms are nonvanishing at 0(0) so that only the qo/T = terms are needed. The 
terms involving g cancel at 0(0) and so only the 0(1) terms are needed (note that 8M± 3> 8M±). Lastly, both the 
0(0) and the 0(1) terms are needed for the terms since 6M± is the same size as 6M±. The first moment is then 
given by 



dn 

4tt ' 



at 



Hq ^ hi 



t| + (i-2/(o)) 5 



Bi 



i - m 



(0) 



1 T\ 



(4.29) 



4 In Ao (eq. [ 1.10 ), the leading order term in e is an even function of u, and the higher order term is an odd function of u. 
Hence, when qg oc ti™ is integrated against Ag to find the moments Mg, the ra = moment is larger than the n = 1, 2 moments 
by a factor of ~ k/m, which are larger than the n — 3,4 moments by a factor of ~ k/m, etc. On the other hand, the leading 
order term in SA± is an odd function of u. The moments of q l against SA±, called <$M™, will then have the n = 0, 1 terms 
larger than n — 2, 3 by a factor of k/m and so on. 
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For the second moment, the B — piece is nonzero at 0(0) so that only lowest order is needed. The B / terms 
must be kept at both 0(0) and 0(1). The result is 



20 



dt 



a/„(k) l = 3 / 1 + 3A 2 



ft,,; Si 



ft, 5, - -%h • B 



(4.30) 



E. Elastic Cross Section 



As discussed before (see §111. A), it is essential to retain the inelasticity in the differential cross section in order to 
derive the correct neutrino flux (as we have done in the previous subsections). Nevertheless, from the energetics point 
of view, v — N scattering is highly elastic near the surface of the proto-neutron star ^ and it is instructive consider 
the " elast ic" sc atteri ng rate, dT/dCl', obta ined by integrating (dT/dk'dSY) over all final neutrino energies k'. Using 
eqs. ( H.ICP and fllllD , together with J4.25Q , we find: 



dT 
dlV 



dk'- 



dT 



dk'dVL' 



M° 



5M°Sl 



B + 5M°_n' B 



( 



Gpcvk\' 



\ 2tt 



1 + 3A 2 + (1 - A 2 )// 



T 



T 



A 1 + A - 4A- n • B + A 1 - A + 4A- ] SI' ■ B 



T 



(4.31) 



The resulting cross-section per particle, (dT/dCl)/n, is similar to eq. (2.4) (recall that polarization P = figB/T for 
nondegenerate nucleons), obta ined by assuming complete elasticity of the scattering process. The difference involves 
the terms with 4AT ' /k in eq. (4.31). These terms appear in the phase space integral of the asymmetric part of the 
cross section, which can be affected by even a small inelasticity. In fact, for low energy neutrinos (with k < 4T), the 
T/k terms dominates the asymmetry in the cross section. ^] 



F. The Anti-Neutrino Cross Section 



The expressions derived in previous sections apply only for neutrinos. For v + N — > v + N , the differential cross 
section of the form eq. ( |4.3| ) still applies, except that one needs to switch the coefficients in front of ft ■ B and fl' ■ B. 
This is due to the crossing symmetry of the tree-level Feynmann diagram. 



V. NEUTRINO ABSORPTION BY NUCLEONS 



In this section we derive an explicit expression the cross-section for neutrino absorption (y e + n — » p + e~) in 
magnetic fields as formulated in §111. B. 

In the regime where the neutrino energy is much smaller than the nucleon rest mass, one might be tempted 
to consider an "elastic" cross-section, obtained by treating the nucleons infinitely massive (i.e., neglecting nucleon 
recoil)Q Further neglecting the effects of electron Landau levels and Fermi blocking, we obtain (e.g., Ref. |53[ ) 



5 In fact, the scattering is elastic (i.e., \qo\ < i) to a good approximation in most regions of the proto-neutron star. The only 
exception is during the first seco nd or so after core collapse, when v e 's are highly degenerate in the stellar core. 

6 Note, however, that eq. (4.31) is valid only for k ^ k m i„ — ^BBim/T) 1 ^ 2 (see the end of §IV.B). This means that the 
maximum asymmetry is ~ \ibB /k m in = (T/m) 1 / 2 . 

7 Previous authors have all neglected the recoil effect. Refs. ^|^3| focus on the rate of neutron decay (including the effect 
of electron Landau levels) and do not address the angular distribution. Ref. [JlTJ discusses the electron contribution to the 
asymmetric emission in e~ + p — > n + i/ e , but the authors did not give an explicit expression. Moreover, the cancellation of 
asymmetric emission and absorption was not considered (see §111. B). 
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cr= -(GWfc) 2 [(l + 3A 2 ) + 2P„A(A + l)fi-B , (5.1) 



7T 



where P n is the neutron polarization, k = ktt is the incident neutrino momentum, and Gf, cy, A are (charge-current) 
weak interaction constants (see Appendix D). This simplified treatment, however, leads to an incomplete expression 
for the asymmetric part (proportional to • B) of the cross-section^. As we show in this section, by incorporating 
the effects of electron Landau levels and small inelasticity, additional asymmetric terms which relate to electron and 
proton polarizations are revealed. These additional terms are important in determining asymmetric flux from the 
proto-neutron star. 

Note that the quantization of electron e nerg y levels can induce oscillatory features in the total absorption cross- 



section as a function of the neutrino energy 60 1 . This effect results purely from the modification of the electron phase 
space due to the Landau levels (similar to the magnetization of an electron gas at low temperatures). This oscillatory 
feature is particularly prominent in the low density regime where only a few electron Landau levels are filled. Our 
focus in this paper is the asymmetric part of the cross-section, which arises from parity violation. We shall therefore 
restrict to the regime where more than a few electron Landau levels are filled. In this regime, we can replace the sum 
over Landau levels by an integral, and obtain an explicit expression for the asymmetric parameter in the cross-section. 

In the calculation presented below, we also include the effect of proton Landau levels. As expected, this introduces 
no new term in the cross-section since many proton Landau levels arc filled for the typical conditions of proto-neutron 
stars. Our result therefore also serves as an explicit demonstration on the validity of using proton plane waves 
in calculating neutrino- nucleon opacities (absorption and scattering; see §VI). It turns out that, after using some 
standard identities involving Landau wavefunctions, the calculation with electron and proton Landau levels is not 
more difficult than the calculation with only electron Landau levels. 



A. The Expression for Absorption Opacity 



The absorption opacity is defined in eq. ( 3.17 ). The energy of a relativistic electron in a magnetic field is given by 

E e = {ml + 2eBN e +pl z ) 1 l 2 , (5.2) 

where N e = 0, 1, 2, ■ ■ ■ is the Landau level index and p e _ z is the electron z-momentum. The other quantum numbers 
specifying the electron states are: er e = ±1, the spin projection along II = p + eA, and R e = 0, 1, 2, • • •, which 
determines the radius of the electron guiding center (see Appendix C). Note that for the ground Landau level, the 
electron spin is opposite to the magnetic field, thus only the spin projection a e o — — sign(p e ^) is allowed. The sum 
over electron states is then 

N e =0<r e =±l R C =0 J - QO 

where c(N e ,a e ) = 1 — <5jv e ,o£cr e ,-cr e o i s zero if both N e = and a e — — cr e ,o and one otherwise. The cutoff R m ax — 
eBA/2ir (the degeneracy of the Landau level) limits the guiding center to lie within the normalization volume V = AL 
(where A is the area). The proton energy is given by 

E p = ^ (^ + £)-W* + f^ (5-4) 

where N p specifies the Landau level, s p = ±1 is the spin projection along the z-axis, fiBp is the proton magnetic 
moment, and p p , z is the proton z-momentum. The summation over states for the proton takes the form 



OO Rn 



rfn ^E E E / L -^r (5 - 5) 

Af=0 R o =0s v =±l J -°° 



8 In traditional laboratory experiments of parity violation (e.g., Ref. [[T2[), the approximation (of neglecting nucleon recoil 
and electron Landau levels) is valid because the temperature T is much smaller than the neutrino energy k. In the supernova 
context, k is comparable to T and this simple approximation leads to significant error. 
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where R p is the quantum number for the proton guiding center, and R ma x is the same as for the electron. Finally, 
the neutron phase space is simply that of a free particle, with energy E n — p^/(2m) — s n /iBnB. 

In Appendix D, we find that the transition rate for absorption (S-matrix squared divided by time) takes on the 
form 

W^ hs) = L- 1 V- 2 (2n) 2 5(E e +E p - k - E n - Q)6( Pe , z + Pp , z - k z - p n>z )\M\ 2 , (5.6) 

where the matrix element, summed over the guiding center quantum numbers R e and R p for the electron and proton 
as well as electron spin er e , can be written as 

R R 

E E E c(N e ,a e )\M\ 2 = ^A^-L^. (5.7) 

R e =Q R p =0 CT e =±l 

Here L uv is the lepton tensor and N^ v is th e nucleon tensor, which takes on precisely the same form as in the zero 
field case. Plugging this back into eq. ( |3.17 ) gives 



AT^CT -00 N P =0 J y ' s n .s p = ±l 

where we have defined a "response function" for absorption 

dpn.z f°° dp PtZ 



-00 27T J_ OQ 2ir 



x (2ir) 2 6(E e +E p -k-E n - Q)6{ Pe , z + Pp , z -k z - Pn , z )f n (l - f P ). (5.9) 

By integrating over the delta functions (see Appendix E), we derive a general expression for the nucleon response 
function: 



ml 1 \ / 1 



ml \e y + 1 J \i + e-y-\ 

where q z = k z -p e . z , qj_ = eB(2N p + 1) - p? i± , q 2 = q\ + qj_, qo =k-E e , 

E n - \i n [i n , Pn,± , [?o + Q ~ q 2 /2m + {^ Bp s p ~ (J, Bn s n )B] 2 ^ Bn s n B 



and 



T T 2mT AT{q 2 z /2m) T 

fj-n - + qo + Q 

z = . 

T 



(5.10) 



(5.11) 



(5.12) 



The above expressions apply for arbitrary values of nucleon degeneracy, recoil energy, and magnetic field. In 
this general case, the nucleon tensor depends on s n and s p , the lepton tensor depends on N e , p e , z , and 
w±_ = |p n ,_L + kx|j and the response function S SnSp depends on s„, s p , N p , N e , p e ,z, and w±. To evaluate eq. (|3.^). 
we are left with two infinite sums and three integrals left to perform. 

To make progress, we shall proceed in the next subsection with an approximate method appropriate to the outer 
layers of the proto-neutron star in which the nucleons are nondegenerate and the recoil energies and nucleon spin 
energies are small in comparison to other energy scales. As discussed before (see §111), only in the outer layers (where 
the neutrino distribution deviates from thermal equilibrium) can asymmetric neutrino flux develops. 



B. Evaluation of the Absorption Opacity: Nondegenerate Nucleon Regime 

For small nucleon spin energies and nucleon recoil, all dependence on s ni s pi and N p can be taken out of the 
exponential in the nucleon response function so that these quantities can easily be summed over. Since S SnS is 
expanded to linear order in (J-bB, it contains only terms linear in s n or s p , but not both s n s p . As a consequence, any 
terms in L^ LV containing s n s p can immediately be dropped, considerably simplifying this expression. We will also drop 
all terms in which will give small corrections to the angle- independent, B = opacity. Lastly, we drop terms in 
L^y which will give zero in the sums over N p (see Appendix D for discussion.) 
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1. Contribution from the N e = State 



Since the electron in the the ground Landau level can only have a spin opposite the magnetic field, the N e = 
term in the opacity expression (eq. |)5.8| ) requires special treatment. As eq. ( |5.8| ) already contains a prefactor £?, we 
can drop all nucleon polarization terms when evaluating the N e = contribution to the asymmetric opacity (This 
cannot be done when summing over all the N e > 1 states since N e = p 2 ± /2eB is summed over a large number of 
states so that the prefactor of B effectively cancels). Since only the nucleon polarization terms contain pieces with 
large coefficients, the N e — state can be evaluated to lowest order in inelasticity. 

In Appendix [e| it was shown that to lowest order in the inelasticity, the nucleon response function for N e = can 
be written 

Ss - sp = exp ij ~ _ u2 ) {Ne = state) ' (5 - 13) 

where u is a dimensionless electron z-momentum defined by 

p eiZ = ±(k + Q)(l + eu), (5.14) 

and 

m J k + Q y ' 

is a small parameter (q z o = k z =p (k + Q)). As this expression for S SnS is independent of N p , we may sum over N p 
in L^ v , with the result (Appendix D): 

oo 

N^L^Ne = 0) = <d(p e ,z)(c 2 v - 4)0 • B. (5.16) 

N p =0 

Since this expression is independent of s n and s p , their sums give a factor of four. The asymmetric opacity from the 
electron ground state is then given by 



2 2tt 2vr v j ' J (2tt) 



,2 



x 4^ exp (f - 1^ - u 2 j e(p e , z )(4 - £)n ■ B 

G 2 F eB^^mT m f^n\,ji _ 2ln „ f°° dp e ^ z 



exp (^) (^ - 4)0 • B /°° ^i(l - /e ) e ^ 2 



2 2^ 2tt \q zfi 

[1 - / e (A + Q)] (4 - 4)H • B, (5.17) 



where we have used eq. ( 4.12) to relate the neutron chemical potential fj, n to its number density n n . In evaluating the 
p e;Z integral in eq. ( |5.17j ), we have approximated f e (E e ) by f e (k + Q) since the first term in an expansion of f e {E e ) 
about this value is odd in u and hence gives vanishing contribution and the second order term in the expansion is 
down by a factor T/m. 



2. Contribution from the N e > 1 States 

For electrons in the excited Landau levels (N e > 1), the relevant matrix element can be written as (see Appendix 
D): 

N^L^(N e >l) = T +T e +T n , (5.18) 

with 
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To 
% 



(4 + 34) 



|A| 



+ 4 e ,AT» ( 1 + 



|A| 



ca(ca + c v )s n ft ■ B - c A (c A - c v )s p Cl ■ B 



(4 - £)n • b, 



7-2 



-l,jV„( w ) 1 



|A| 



7^ h(i + p ' 



|A| 



(5.19) 
(5.20) 

(5.21) 



where |A| = (p 2 ez + 2eBN e ) 1 / 2 , oj = w 2 , I (2eB) = (p n ,X + kj_) 2 /(2ei3), and the function In c ,n p (The shape of the 
Landau wave function) is defined by eq. (|C4j). Note that in these expressions, we have dropped all terms that will 
give zero contribution to the opacity (such as those terms invo lving s n s p ). 

The N e > 1 contribution to the absorption opacity (eq. |5.8|) can be written as 



K (abs) (iVe > j j = ^(abs) + K (M)( e)Ne > !) + K (abs) ( „ p); 

where each piece corresponds to a different part of L^(N e > 1): 

(abs) _G 2 F eB ^ r°° dp e , z ^ r d 2 p n ^ ^ 

)\r 1 w — OO \r n « \ / ^ ^ 



jV e = l 

k ..m m . 3: g f; /~ f^ (1 _ A) £ /^i E s ,„„r„ 



N„=0' 
00 ,, , 2 



(5.22) 

(5.23) 
(5.24) 
(5.25) 



We are going to evaluate the sum by replacing it with an integral. Such a procedure effectively eliminates 

any possible oscillatory behavior of the opacity as a function of energy (see the beginning of §V) , but is valid when 
more than a few electron Landau levels are filled. For infinite nucleon mass, energy conservation r equire s k + Q = 
(p 2 z + 2eBN e ) 1 / 2 (neglecting m e ). For given N e and p e ,z, the nucleon recoil energy is of order \q z \ yjT/m. Thus it is 
natural to define a dimensionless recoil energy u via 



2eBN e = E\(\ + eu), 



where 



{k + Qf-p 



2 

e,2 ; 



8T\ 1/2 \q z \(k + Q) 



The nucleon response function (eq. 5.1C ]) can be expanded for small eu, with the result (Appendix E): 

rn 



S SnSp = 1 — r exp(-y )(l - Sy), 

\q*\ 



where 



and 



, s ( Pn,± 

exp(- yo ) ~ exp - - _ 



eu 6 E 2 ± 
2{k + Q) 2 



2u(k + Q)q 2 



emE 2 ^ 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



Sy 



l^BnSnB 

2T 
HBpSpB 
2T 



tumE\ 



euE' 2 



<& (k + Q)q 2 V 4(fc + Q) 2 

q 2 , tumE 2 , ( tuE 2 , 

1 + — H — 1 — 

^ q 2 + (k + Q)q 2 \ 4(fc + Q) 



(5.30) 



Similarly, for small eu, the electron Fermi blocking factor can be expanded to first order in e as 
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1 + ^Tg)^ + « 



(5.31) 



Now consider K^ h ^ in eq. ( 5.23 ). Since only S SnSp depends on spin, the spin sums J2s n s effectively set Sy = 0. 
The factor exp(— yo) can be evaluated at lowest order. The sum J2n can ^ e calculated using the summation rule 
for Ijf S (Appendix D). Replacing J2 N by J dN e = J eEj_du/(2eB) and integrating over p n> _\_ (see Appendix F), we 
arrive at 



(abs) 



Gl eB 



mT 



(k+Q) 

-(k+Q) 2?r \Qz\ J-ao 

V*(hj-n\ r(k+Q) 



2 2tt v v a > 2tt \m J 2eB 

^(k + Q) 2 n n (c 2 v + 3c 2 A )[l - f e (k + Q)], 

7T 



dpe.z Tn f°° eE^du 
2eB 

dp, 



~(k+Q) 



2tt 



exp(- U 2 )(l-/ e ) 

(We) 



(5.32) 



which is exactly the usual B = opacity. 

Next consider K^ bs \e, N e > 1), the "electron contribution" from the N e > 1 states to the opacity (eq. [ 5.24 1). 
We may evaluate all quantities to lowest order in the inelasticity. The spin sums and nucleon response function are 
evaluated as before. Performing all integrals but p eiZ gives 



r(k+Q) 

^ abs )(e,A e > l)oc / dp e . zPe 
J -(k+Q) 



0. 



(5.33) 



So the electron contribution from the higher Landau levels is zero to lowest order in the inelasticity. The next order 
correction scales as T/m, and c an be neglected. 

Finally, n^ hs \np) (eq. ||5.25| ) gives the contribution of nucleon polarizations to the opacity. Performing the spin 
sums, the integral over p n ,±, and using the results of Appendix D for the sum over N p yields 



^e^fl-fiT ^™ f>-/ e )exp( 
27r J -co 2tt \q z \ ^ 



HBnB ^BpB 

- < T c A (c A + c v ) ^—ca(ca-c v ) + 



VBnB ^BpB 

T CA(CA + C V )-\ 7^—Ca{Ca - C V ) 



e u m 



E\ 2eBN e + k 2 ± +eBp e>z /\A\ 



4(fc + Qfq 2 z 

€umE\ E 2 u i mE 4 L 



(k + Q)q 2 2{k + Qfq 2 



2u 



2 q 2 + 2eBN e + k\ + eBp CtZ /\A\ 

q 2 



(5.34) 



Again changing the sum over N e into an integral over u, expanding / e , and then performing the trivial p e z integral 
gives the final result 



G 2 , 



«l»MJ(„p) = r±l( fc + Qf nn n . B [1 - f e (k + Q)] i 2^^ca(ca + c v ) 

IT 1 



,B 



T 



{k + Q) 



(fc + g) . 

1 + f fe{k + Q) 



2 c A (c A + c v ) + 2 £ c A [c A - c v ) 



(5.35) 



3. Result 



We can summarize our result for the absorption opacity in the following transparent formula: 

K (abs) = K (ab S ) ( 1 + eabsf2 . B ) ; 

where k^° s ^ is the B = opacity as given by eq. ( 5.32j ) , and the asymmetry parameter e a b s is given by 



(5.36) 
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e a bs = £abs(e) + e abs (np), 



(5.37) 



with 



e a bs(e) 



r 2 
\2 A 



e a bs(np) 



1 eB 



2 (k + Q) 2 c 2 v + 3c 2 A 

Ca(ca + Cv) M-Bn-B 



34 



(fc + Q) 



1 



(*• 



T 
Q) 



/e(fc + Q) 



C A (CA + Cy) HBnB Ca(ca - C V ) fJ-BpB 



3c 2 A 



T 



3ci 



T 



(5.38) 



(5.39) 



Comparing this result with eq. ( |5.lD , which was obtained from a simplified calculation assuming infinite nucleon 
mass and neglecting Landau levels for electrons, we see that the simplified calculatio n g ave an incorrect result for 
the asymmetr y pa r amet er. Only a neutron polarization term was included in eq. (5.1). The correct expression 



for e a bs (eqs. |ft.37| -| 5.3S |) contains an electron contribution (arising from from the ground-state Landau level) and 



additional contribution from both neutron and proton polarizations (which arises from our more careful treatment of 
inelasticity). These news terms dominate the asymmetric parameter for neutrino energy k/T< a few. It is interesting 
to note that the particles in the final state of this reaction (p and e~) contribute to the asymmetry parameter. 
Previous investigators found a contribution only from the initial state particles. 



C. Moments of the Absorption/Emission Rate 



With the absorption opacity given in previous subsection, it is straightforward to calcul ate t he moments of the 
abso rptio n/emission rate in t he Bo ltzmann equation. The abso rption rate is given by eq. (3.15) with 5f u given by 
eq. ( 3.21 ) . Using the opacity ( 5.36 ) and the integrals in eq.( 4.19 ), we obtain the zeroth, first and second moments: 



47T 



-p. . 

4tt y 



dfuM 



at 
a/,(k) 



dt 
%(k) 



dt 



abs 



abs 



abs 



^*(abs] 



*(abs) 



* (abs 



+ e a bsh 



•B , 



1-2 



hi + Y^fabs ( Bihj + Bjhi - ~<S;,h ■ B 



(5.40) 
(5.41) 
(5.42) 



where 



*(abs) _ (abs) 
K Q — K Q 



1 + exp 



T 



(5.43) 



(see eq. fig]). 



D. Absorption Opacity for v e 

So far we have been concerned with the opacity for i> e + n — ► p + e. The opaci ty for D e + p — > n + e + can been 
obtained by a similar calculation. The result can be found from the equations in §VB3 by replacing Q by —Q, n n 
with n p , the electron distribution function by the positron distribution function (which is rather small), and nsn with 
—fiBp- Thus, the absorption opacity for v e is given by 

s (abs) = -(abs) ^ + - absn . ^ (5 44) 

where Kg abs ^ is the B = opacity as given by 



H„ 



(abs) 



G 



2 



„. \ n "W '"V 

TT 



(k - Q) 2 n p {c 2 v + 3ci)[l - f e+ (k - Q)]0(k - Q), (5.45) 



and the asymmetry parameter e a b s is given by 
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^abs — ^abs(C 



e a bs(np), 



(5.46) 



with 



e a bs(e + ) = 



1 eB 



2 (k - QY 4 



3ci 



e a bs(«p) 



r 



T 



(fc-Q) 



(fc-Q) 

T 



fe+(k-Q) 



t C A (c A - Cy) [i Bp B n C A (c A + Cy) /XB n .B 



3ci 



r 



34 



r 



(5.47) 



(5.48) 



The theta function in eg. ( |5.45 ) comes from the fact that the reaction is not energetically allowed unless k > Q since 
the proton is lighter than the neutron. 



VI. THE MOMENT EQUATIONS OF NEUTRINO TRANSPORT 



In the last two sections (§IV and §V) we have carried out detailed calculations of neutrino scattering and absorption 
in magnetic fields. Explicit expressions have been obtained in the nondegenerate nucleon regime, which is appropriate 
for the outer layer the of proto-neutron star where asymmetric neutrin o flu x (drift flux) is expected. We now use these 
results to derive the moments of the Boltzmann transport equation (3.1). We focus on v e below, although similar 
results can also be obtained for other neutrino species. 



The zeroth moment is obtained by integrating eq. (3.1) over dfl. Combining eqs. J426| ) and ( p5| ) we find 



(0) 



at 



r dn 


W„(k)" 




r dn 


ra/„(k)i 


1 4^ 


at 


+ 

sc 


J 4^ 


at 



1 1 ■ - 



k 



T 



dg 
dk 



B - K, 



^(abs 



g + e abs h • 



(6.1) 



The zeroth moment equation governs the energy exchange between matter and neutrinos. The B = part of the 
scattering opacity can be ignored in the zeroth moment equation since it is suppressed by a factor of k/m (i.e., only 
the inelastic part of the scattering contributes to matter-neutrino energy exchange; recall that we have not included 
inelastic electron-neutrino scattering which can be a much larger effect). The term — K^ ahs ^g represents the usual 
neutrino emission and absorption. It is of interest to note that the asymmetric parts of the scattering and absorption 
introduce new terms to the zeroth moment equation. The importance of these terms will depend strongly on the field 
strength and optical depth. In a steady state (neglecting the time derivative term) we have 



Vh 



*(abs) 
l ■< 



*(abs), -p, (sc) 



T dh 4 T h 

dk k 



B. 



(6.2) 



The first moment equation is obtained by multiplying eq. (3.1) by £li and then integrating over d£l. Combining 
eqs. ([1.29|) and (|5.41[) we find 



dt 



dxi 



-Kq ^ hi 



3 dxj 

(sc) 
^sc^o 



dn { 

47r' 



0/,(k) 



dt 



t|| + (1-2/(°)) 9 



Bi 



^absgBi + -e a bs!ijBj 



4tt 



dt 



abs 



1 - 2/W + 1 - 4^ 1 J,,B 



(6.3) 
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The time derivative term can almost always be dropped (This corresponds to a rapid redistribution of matter tem- 
perature, the timescale of which is of order the mean free path divided by c, much smaller than neutrino diffusion 



time of the star; see Ref. j73|). Defining the total B 

1 d(ti 0) +g) : 



opacity, Kg tot ^ 
din 



*(abs 



we have 



3 k 



(tot) 



dxi 



(tot) Q X 



3 k 



(tot) 



>c) 



3k l 



*(abs) 



Eabsff r & 



3 k 



.(tot) 



2k, 



= (abs) 



^abs 



(6.4) 



Clearly, in addi tion to the usual diffusive flux (the first line of eq. [3.4]), there is also a drift flux (the second and third 
lines of eq. [6.4]) which depends on the direction of the magnetic field. This asymmetric drift flux is a unique feature 
of parity violation in weak interactions. Equation ( |6.4[ ) explicitly shows that the drift flux is nonzero only when the 
neutrino distribution deviates from thermal equilibrium, as expected from general consideration of detailed balance 
(§111. C; see also §11 for discussion). The spherical deviation, g, always gives rise to a drift flux along B. However, the 
drift flux from Jy is along the direction of the vector IijBj, which does not have to be directed along the magnetic 
field at all points in the star. For cylindrical symmetry, however, one would expect that the net flux produced by the 
IijBj term would average to the B direction. 

Finally, the second moment equation can be obtained by multiplying eq. (3.1) by Vij = (3Qjfij — <5jj)/2 and then 
integrating over d£l. Combining eqs. (4.30) and (5.42) and ignoring the time derivative, we find 
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dt 
-2/i° 



1 
A 



47T ' 

. T 



3/,(k) 



Of 



hiBj 



abs 

hjBi 



*(abs) 



Ii 



3 

To 



£abs ( Bihj + Bjht - -<5 y h • B 



(6.5) 



The above equations apply to v e . Similar equations can be derived for other species of neutrinos. Note that 
since &V T ) and £v T ) are always created in pairs inside the proto-neutron star, they have the same energy density 
distribution. Because of the crossing symmetry (see §IV.F), the drift flux of v^M exactly cancels the drift flux of 



VII. DISCUSSION 



In this paper we have presented a detailed study of neutrino-nucleon scattering and absorption in strong magnetic 
fields. Specifically, we focused on the effect of parity violation in weak interactions which can induce asymmetric neu- 
trino transport in the proto-neutron star. Starting from the weak interaction Hamiltonian, we found the macroscopic 
moment equations of neutrino transport. Explicit results applicable to the outer region of a proto-neutron star are 
given in eq.s (6.2), (3.4) and (6.5). Despite the fact that the neutrino cross-sections are asymmetric with respect to 



the magnetic field throughout the star, asymmetric neutrino flux can be generated only in the outer region of the 
proto-neutron star where the neutrino distribution deviates from thermal equilibrium. 

Previous studies based on simplified treatments (see §11) have led to misleading results. We have tried to clarify many 
of the subtleties in deriving the correct expressions. The main technical complication lies in the proper treatment 
of the inelasticity of neutrino-nucleon scattering/absorption: although these processes are highly elastic from the 
energetics point of view, it is essential to include the small inelastic effect in order to obtain the correct expression 
for the asymmetric neutrino flux. In addition, it is necessary to use Landau wave functions for the electron since the 
quantum mechanical ground state of the electron gives the dominant contribution to the asymmetry for low energy 
electron neutrinos. To obtain simple formulas for the respective opacities, we devoped a method to expand phase 
space integrals for both small magnetic field strengths (and correspondingly small spin energies) and also for small 
inelasticity. This method has general applicability for computing the effects of nucleon recoil in phase space integrals 
in powers of T/m. 

To quantitatively determi ne th e asy mmet ry in neutrino emission from a magnetized proto-neutron star one has to 
solve the moment equations (|6.2|), (EOI) and (tLa) in the outer layer of the star. This is beyond the scope of this paper. 
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Here we shall be contented with an order-of-magnitude estimate. First note that the net drift flux associated with 
VfiM and Sv T ) is zero, and we only need to consider v e and v e . The key quantities that determine the asy mmetr ic 



flux are the dimensionless asymmetry parameters e sc and e a b s . For neutrino- nucleon scattering, we have (eq. | 4.28 |) 

e sc ~ 0.006 B lb T-\ (7.1) 

where -B15 is the field strength in units of 10 15 G, and T is the temperature in MeV. The asymmetry parameter for 
neutrino absorption has contributions from electron and nucleons: 

£abs (e)~0.6 5 15 £; 2 , e abs (n.p)~ 0.006 SigT" 1 [1 + 0(T/E U )}, (7.2) 

where E v is the neutrino energy in MeV. Thus for high energy neutrinos the asymmetry is dominated by e sc and 
£abs(fip) (arising from nucleon polarization, ~ pbB/T), while for lower energy neutrinos it is dominated by e a bs(e) 
(arising from electrons in the ground Landau level). The electron neutrinos decouple from matter near the neutri- 
nosphere, where typical density and temperature are p ~ 10 12 g cm~ 3 , and T ~ 3 MeV. For a mean v e energy of 
10 MeV, e a b s is greater than e sc . The asymmetry in the v e , v e flux is approximately given by the ratio of the drift flux 
and the diffusive flux, of order e a b s [KQ abs V K o t0t ^]- Averaging over all neutrino species, we find the total asymmetry 
in neutrino flux a ~ 0.2e a b s - To generate a kick of a few hundreds per second would require a dipole field of order 
10 15 - 10 16 G. 

Since the asymmetric neutrino flux depends crucially on the deviation of the neutrino distribution function from 
thermal equilibrium, it is of interest to consider how the function g = f v — f„ (or Jy) scales with the depth of the 
star measured from the surface. Without a magnetic field, we expect g to decrease exponentially toward zero below 
the decoupling layer (which is close to the neutrinosphere for v e and v e ). In the presence of asym metric absorption 
and scattering opacities, this scaling may be modified. Inspecting the zeroth moment equation ( |6.2|) , we may conclude 
g ~ eh ■ B in the deep interior of the star (Recall that in radiative equilibrium one always have V • jh.dk = 0). This 
effect may increase our estimate for the asymmetric flux. We hope to address some of these issues in a future paper. 
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APPENDIX A: THE MATRIX ELEMENT FOR v-N SCATTERING 

In the usual case in which the spin projection of the particles is not "measured" , one can sum the matrix element 
over the final spins and average over the initial spins. However, for spin 1/2 particles in a magnetic field, the initial and 
final state spin dependence in the matrix element cannot be immediately summed over since the nucleon distribution 
function and the energy conservation delta function both have spin dependence of the form —sfigB/T. To obtain the 
differential cross section one needs to calculate the matrix element \M ss >(n, fi')| 2 for initial (final) nucleon spin s (s') 
and initial (final) neutrino direction Jl (fi')- We neglect the effect of Landau levels of proton in the v — N scattering 
cross-section (but see Appendix C and D). 

The low energy effective Hamiltonian density for neutral current scattering of a spin-1/2 fermion with a neutrino 
is given by (see, e.g., Refs. |Io| , |7^] ) 

Hint = (CV - CA15) ^N^u'f (1 " 75) + h.C, (Al) 

where neutral current vector and axial coupling constants are |lC| ] given by^| cy — —1/2 and ca = —1.23/2 for 
!/ + n->i/ + nandcy = 1/2 - 2 sin 2 9 W = 0.035 and c A = 1.23/2 for v+p^ v+p. Here G F = 1.166 x lO^GeV" 2 
is the universal Fermi constant and sin 2 9w = 0.2325 (9w is the Weinberg angle. 



9 Raffelt and Seckel |74|,(75| considered the isoscalar contributions to the scattering amplitude as well as the usual isospin pieces, 
and suggested cy — —1/2 and ca = —1.15/2 for v + n — > v + n and cy = 1/2 — 2 sin 2 9w and ca = 1.37/2 for v + p — > v + p. 
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The (nonrelativistic) nucleon wavefunction with four momentum P = (m + E, p) ~ (m, 0) and spin four-vector 
S ~ sB is given by 



= V- 1/2 U N e lp - x -' lEt , 



(A2) 



(where V is the normalization volume and Un is the 4-spinor) , while the neutrino wavefunction with four momentum 
K = (k, Ml) is 



^ = v --l/2 [ ^ e tk.x-tfct < 



(A3) 



For the antineutrino, replace U v exp (ik • x — ikt) with exp (— ik • x + ikt). 

Plugging the wavefunctions into eq. QA1| ), the transition rate W (S-matrix squared divided by time) can be written 



where 



w ({ P , s, si} , { P ', s ', n'» = (2^) 4 5 4 (p + K-pi- k') \m ss , (n, n')l 2 



(A4) 



(A5) 



V - tf'„7„(l ~ H>)U V U V ~( V (\ - 15 )U' V = — Tr [JC' 7 „(1 - 75 ) £>(1 - 7 s)] 



—K' a K p X a ^, 



Xarfu = -Tr [ 7q7 ^ 7/ 3 7! /(1 - 7 5)] = 3^5/^ _ 9<*p9nu + 9av9Pv, - itanpu, 



(A6) 
(A7) 



N^(s,s') = U' Nl »{c v - c Al5 )U N U Nl »(c v - c Al5 )U' N , 



(A8) 



and we use the sign conventions goo — +1) 9ij — an d £0123 = +1- The nucleon piece can be evaluated using the 

spin projection operator |7(J (1 + 7s73s)/2 and the energy projection operator (1 + 7 o)/2 so that 



N^ u (s, s') = -Tr [(1 + 7573s') (1 + 70) l^{cv - caJs) (1 + 7 5 7 3s) (1 + 70) Y(c v - ca 7 s)] 
Explicit computation of each component gives 

AT° = ~4(1 + SS ') 

N °i = N io = _l A ( s + s > )§ . 3 



(A9) 



Nij = 2 C a - ss ') + 2ss'5 i3 5 j3 + ie 0ij3 (s' - s)] 



(A10) 



where e^ u \ a is the completely antisymmetric tensor with 60123 = +1. The remaining traces can be evaluated by 
standard methods 76 7% with the result 



|M„/(n, n')\ 2 = 1 -g 2 f c 2 v { (1 + 3A 2 ) + (1 - a 2 ) n ■ n' 



-2A(A + l)(sH + s'ft') ■ B - 2A(A - l)(sSl' + s'Sl) ■ B 

(1 - a 2 ) (1 + n • n') + A\ 2 n ■ bo' • b j 



(All) 



where we have defined A = c^/cy. 

Time- reversal invariance can be explicitly checked for the matrix element in eq. (All), or equivalently the S-matrix 
in eq. {A4), by simultaneously exchanging all initial and fina l sta te labels. 

For antincutrinos, one would just switch and f2' in eq. (All). 
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APPENDIX B: NUCLEON RESPONSE FUNCTION FOR SCATTERING 



Following the procedure outlined in [Q, we first use d 3 p' to integrate over S 3 (p + q — p') and then integrate over 
the azimuthal angle for p, with the result 



f f°° f 1 

Sss'(qo, q) = ^ J dpp 2 J ^ d^S(q + E- E')f N (E) [1 - f N (E')} , 



(Bl) 



where fi = p • q/pq and E' = —fisBs' + (p + q) 2 /2m = —[IbBs' + (p 2 + q 2 + 2pqfi)/2m. Care must now be taken to 
correctly integrate over the energy-conservation delta function with the coordinate /x; the integral is only nonzero if 
the argument of the delta function is zero for pi = ^ € (— 1, !)■ We find 



Mo = 



so that Hq < 1 for 



2 \ 2 
5 >Pmin 



qa - q 2 /2m - fj, B B(s - s') 
pq/m 



q - q 2 /2m - fi B B(s - s') 



-i 2 



q/m 



Changing variables from p to E = —sfigB + p 2 /2m in the remaining integral gives 

,2 i-oo 



m 2 Z 100 

S ss <(g , q) = ^— dEf N {E) [1 - M£ + q ) 



where 



'M2 



2m 



= -[IbBs 



[ g0 - g 2 /(2m) - [i B B{s - s')] 
4(g 2 /2m) 



(B2) 



(B3) 



(B4) 



(B5) 



is the minimum energy allowed for the initial state nucleon in order for energy and momentum conservation to 
be satisfied given go,g, s, and s'. The last integral can be done by first defining the dimcnsionless variables x — 
(E - ijln)/T, x mhl = (E min - hn)/T, and z = q /T, yielding 



S ss '{qo,q) 



m 2 T 



dx 



27r <7 Jx„ 
_ m 2 T 1 
27rg 1 — er 



1 



1 



■In 



1/ V 1 + e ~ x ~ 
1 + cxp(-x min ) 



1 + exp(-X m in ~ z) 



(B6) 



This expression agrees with |l0[ for B = keeping in mind that our definition of S is a factor of two smaller than 
theirs. 

Expanding x m j n to linear order in B we find 



■^min — ^0 "T" 3% 

(go - q 2 /2m) 2 _ pi N 
T 

2mq 



x 



Sx 



4T(g 2 /2w) 
-H B B 



2T 



1 



s+ 1 



2mg \ , 



(B7) 



For fa « 1, S ss > can be written as a sum of So, the zero field value, and SS SS ', the correction due to the magnetic 
field, i.e., 

S ss >{q ,q) = S (q ,q) + SS ss >(q Q ,q) 

m 2 T 1 1 + cxp(-a;o) 



So(qo,q) 



2nq 1 — e 
m 2 T 



■In 



1 + exp(— xo — z) 
Sx 



2nq [exp(x ) + 1] [1 + exp(— x — z)] 



(B8) 
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Note tha t the asymmetry in the the coefficients in eq. ( |4.8| ) is entirely due to the 2mqo/q 2 terms, which first appear 
in eq. (BT) as a consequence of the energy and momentum conservation delta function. Had one initially set m — > oo 
in S(E + qo — E'), these terms would not have appeared. 

In the limit of nondegenerate nucleons with (j,n/T -C —1 and cxp(un/T) = (2 1 ' 2 -K 3 ' 2 n)/ (m s ' 2 T 3 ' 2 ), we find 



In 



1 + exp(— xq) 



1 + exp(— z - z) 
'. exp (—xq) — exp 



1 



[exp(x ) + 1] [1 + exp(-£ - z)] 
fJ-N (qo - q 2 /2m) 2 ' 
T 



AT(q 2 /2m) 



(B9) 



which takes the form of a Gaussian in k 1 . The center of the Gaussian is located at k' ~ k and has a width of order 
(T/m) 1 / 2 fc due to the recoil motion of the nucleons. We can further simplify this expression by defining the small 
quantity e and the dimensionless variable u by 



e = 


[4(1 


-H')T- 


1/2 








m 


) 




so that k' = k(l + eu). Then the recoil momentum is 






q 2 = 


2fc 2 (l- M ')(l + 


eu) 


and 










(q - q 2 /2mf ^ ^ 


2 2 
U 


emu 




AT(q 2 /2m) 4T[2fc 2 (l - 


M')/2m](l + 


eu) 


The criteria for the expansion of the B = 


- part are 


T < 


m 



k' 



l-eu + 2 



k(l-u') 



(BIO) 



(Bll) 



(B12) 



require the additional assumptions that both \isEjT -C 1 and (ubB /T)(yTm/k) <C 1. To first order in e we then 
find 



So(qo,q) 



ek 



■ exp(— u 2 ) 



1 % 
1 - -eu + eu 6 - 2 



k(l-u') 



and SS SS ' = —SqSx involves the quantity 

2mqo 



k(l-u') 



[u — eu 



(B13) 



(B14) 



The range of the variable u is from u min — — l/e<C— 1 to oo. Since e«l we may extend the lower limit to — oo with 
only exponentially small error. 



APPENDIX C: WAVEFUNCTIONS FOR THE ABSORPTION OPACITY 



The electron wavefunctions in cylindrical coordinates (p, (j>, z) are given in |7q |. In the standard representation for 
the Dirac matrices |f76f , these wavefunctions are written as 



* e = i" 1/2 e ip -*- <B '*^ B (p,0) , E e = ( m 2 e + 2eBN e +pl z ) 1 / 2 (CI) 
where p e z is the z-momentum, and N e = 0, 1, 2, • • • is the Landau level index. The 4-spinor U e is given by 

1/2 



Ue(p,4>) 



1 



2ttA 2 



( C 1 I Nc ^ Rc (t)e-^\ 

iC 2 lN e ,R e (t) 

C 3 I Ne - llRs (t)e-^ 
\iC 4 I NetRe (t) J 



(C2) 



with 
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Ci = a + A + 
C 2 = a e a-A^ 
C 3 = cr e a + A- 



a± 
A± 



(C3) 



!/i = (5) 12 e- i/2 i (n - r)/2 4"- r) («) , / dll'lil, I. (C4) 



7f! 



Here A = (ei?) -1 / 2 is the cyclotron radius, i = p 2 /2X 2 (not to be confused with the symbol for time), and V = LA is 
the normalization volume (with length L along the z-axis and area A in the x-y plane; one particle in volume V"). The 
operators which have been simultaneously diagonalized, and their corresponding eigenvalues, are perpendicular energy 
Ej_ = (eB /m e )(N e + 1/2), z-momentum p ejZ , the radius of the guiding center R gc — X(2R e + l) 1 / 2 (R e — 0, 1, 2, • • •) 
and the "longitudinal" spin polarization operator A = E • II = cr e (p 2 z + 2eBN e ) 1 / 2 , where II = p + eA and cr e = ±1. 
Note that A reduces to helicity when B = 0. The coefficients Ci are found by requiring the spinors to satisfy the the 
Dirac equation, the normalization condition, and by diagonalizing A. The functions I nr {t) are the same as in [ [78[ , 
except that we have written them in terms of the Laguerre polynomials as defined in |79| . 

As the properties of the ground state are crucial to our results, we now describe them in some detail. The lowest 
energy state has N e = 0. As the functions I-x,R e are not well defined, only the spin projection a e ^ — — sign(p e ,z) 
is allowed for N e — 0; hence the coefficient of I-i,R e will be zero. To interpret this restriction, remember that a e is 
the projection of the spin on the vector it, not the magnetic field. However, by expressing 7r in terms of raising and 
lowering operators one may easily show that (E x ir x + E 2/ 7r 2/ )^ r (AT e = 0) = 0, and hence A x S(N e = 0) = S^7r z vI/(A e = 
0) = Yj z p z ^(N e = 0) = — p e>z i if(N e — 0). Hence the N e = state with a e = a e fl has spin opposite to the magnetic 
field, as expected of the electron ground state. Since only one value of a e is allowed for N e = 0, this state is expected 
to have a different form for the matrix element than the N e > 1 states. 



We assume that the protons are non-relativistic in which case the wavefunctions can be written as 76 

* p = L-^e^'-^Upip, (/>), (C5) 

where 



( S Sp . +1 I Rp . Np (t) 

8 Sp ,-llR p ,N p (t) 



Vo 



(C6) 



The energy for the proton is E p = (eB/m)(N p + 1/2) — [Ab p Bs p + p p z /2m. When there is no subscript on the mass, 
we have approximated m n ~ m p = m. The spin projection along the magnetic field is s p = ±1. The non-relativistic 
spinors differ from the relativistic spinors in that they are completely decoupled from each other. Furthermore, the 
different components of the spinors correspond to different energies due to the anomalous magnetic moment (Recall 
that in the electron case, one could arrange the spinors so that the energy depends only on N e ). This complicates the 
calculation of the matrix clement since the proton distribution function will now depend on spin (through the energy) 
and the matrix element cannot be directly summed over. Lastly, note that the order of N and R was switched going 
from the electron to proton case, since the electron has angular momentum L z (e) = N e — R e and the proton has 
angular momentum L z (p) = R p — N p due to the sign of the charge. The ground state of the proton is the state with 
N p — and s p = +1. Both s p = ±1 are allowed for N p — 0. 

The neutron and neutrino wavefunctions are the same as those used for the scattering calculation (Appendix A). 

APPENDIX D: THE MATRIX ELEMENT FOR ABSORPTION 

The S-matrix is given by 

S fi = -i J d 4 xH mt (Dl) 

where the weak interaction, low energy effective Hamiltonian is (Refs. []l0| , [74| ) 

G 

Hint - ^/§*p7m (cv ~ c Alb ) *„W e7 ^ (1 - 75 ) + h.c. (D2) 
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where the neutral current coupling constants for the absorption process are |74] cy — 1.00 and ca — 1-26 (we shall use 
the same notation for these coupling constants for absorption and scattering, eve n though their values are different). 
Plugging the wavefunctions from Appendix C and the Hamiltonian in eq. (|D2]) we find 



G 

S f i = -i-^L^V-^nSiEe + E p -k- E n - Q)2ir5{ Pe , z + Pp , z -k z - Pn , z ) 
v2 



x / dpp 
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U P (P, 0)7/i ( c v - CAlb) U n U e (p 1 0)7 M (1 - 75) U„ 



(D3) 



where wi = ( Pn ,x + k x )e x + (p n . y + k y )e y , and xj_ = xe x + ye y . The transition rate (the square of Sfi divided by 
time) can be written in the form 



W- f 

where we have defined 

r 2 

\M\ 2 = ^ 



(abs) _ I Sfj 1 2 _ J 



= L- 1 V- 2 (2nYS(E e +E p -k-E n - Q)5{p e>z + Pp , z - k z - Pn , z )\M\ 



(D4) 



2tt 



dpp / dxt>e™^U p {p, 0) 7/J (c v - c Al5 ) U n U e (p, </>)-f (1 - 75) U v 



(D5) 



The integrals over p and </> can be accomplished using eqs. (4.6) and (4.7) in Ref. jf8|, which in our notation give 

1-2-K 

2tt 

where J/v(z) is the n-th Bessel function, and 



(D6) 



dpp 



^lN 1 RAt)lR,NMJ(N 1 -B, 1 )-(N 2 -B, 2 )(^2t\w 1 _) = {-l) N ^ R H NlN2 {\ 2 W 2 j2)I RlR2 {\ 2 W 2 j2). (D7) 



After performing these two integrals, the matrix element \M\ 2 can be written as 



\M\ 2 = ^fl ReRp (X 2 w 2 j2) C/ p7 m (cv ~ c Al5 ) U n U e r (1 - 7s) U v 



where we have defined 



Uu 



( S s p ,+i 


Vo 



and U e 



/C 1 / JVc _ 1 . A r p (A 2 W 2 L /2) 

iC 2 I Nc . Np (\ 2 wi/2) 
C 3 I Ne . ltNp {X 2 w 2 j2) 
\iC A I Ne , Np {X 2 w 2 ± /2) 



(D8) 



(D9) 



At this point, the summation rule for the I ns (eq. [2.24] in Ref. ifTq] ) may be used to sum over the guiding center 
coordinates with the result 



Rmax Rn 



E E'U=E 1 = ^ 



R c= i? B =0 



-R„=0 



eB_ 
2?r' 



(D10) 



Since only the matrix element depends on a e , it may be directly summed over. As a result, the matrix element may 
be put into the form 



Rmax Rmaoi 



Gl .eB 



E E E c(N e ,a e )\M\ 2 = ^fA^ Lflv N^ 



i? c =0i?. B =0<T c =±l 



(Dll) 



where the lepton piece is 



L ^ = E C (^ V e,C r e)?7 e 7 AI (1 - 75) U V U V ^ V (1 - 75) U e 



<T e =±l 



— 53 c(N e , a e )Tr \u e U elfl (1 - 75) (1 - 75) 



a.=±\ 



(D12) 
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and the nucleon piece is 



N"" = U 



r p7 M (cv - cats) U n Unl v (cy - ca7s) 



(D13) 



A moment of inspection shows that the form of the nucleon tensor is exactly the same as eq. (A1C) for the scattering 
problem if we replace the initial nucleon with the neutron and the final nucleon with the proton. The dependence on 
the shape of the proton wavefunctions is contained entirely in the In c ,n functions. 

The lepton tensor is more complicated. In particular, notice that the N e — electron Landau level is the only state 
for which there is only one polarization. Consequently, the matrix element will have quite a different structure for 
the electron ground Landau level and will contain the important electron contribution to the parity violation effect. 
The answers for N e = and N e > 1 will be given separately. 

In the N e = case, by representing U e U e in terms of gamma matrices, performing the traces, and then summing 
against N 1 *", we find 

L^N^(N e = 0) = &{Pe, z )ll Np (w 2 j2eB) [<$ + 3c 2 A + {c 2 v - c 2 A )Sl ■ B 
+2c A (c A + c v )(s p + s n n ■ B) - 2c A (c A - c v )(s n + Spil ■ B) 



+s 



nSp {4 -C 2 A + ( C 2 y + 3c\)n ■ B} 



e(Pe, z )ll Np «/2eB)(c 2 v 



c 2 A )n-B, 



(D14) 



where we have taken the relativistic limit (m e — » 0) for the electrons. For simplicity, several terms in this expression 
which do not affect the final result for the opacity have been discarded. First, we have thrown away terms which 
will only give corrections to the angle-independent piece of the opacity. Second, any terms with s n s p have been 
dropped since they will always give zero in the sum over s n and s p (since we are expanding the nucleon response 
function to linear order in the spin energies). Finally, we have dropped terms like s„f2 • B and s p £l ■ B since the 
N e = contribution is already proportional to B; these terms will yield an additional factor of B from the nucleon 
polarization which will be much smaller. 

We shall need L^N^ (N e — 0) summed over all N p . Using the summation rule for Ins'- 



s=a 



we find 



00 

E 

N„=0 



L^N^iNe = o) = e(p e , z )(4 - 4)n ■ b, 



(D15) 



(D16) 



where & is the step function. 

For the N e > 1 case, all terms containing a factor s n s p , which gives zero in the sums over s n and s p , and corrections 
to the angle independent opacity will be dropped. Furthermore, we drop terms with no spin dependence which 
are proportional to In c .n -T/v c -i,jv ■ Since these terms have no spin dependence, they cannot couple to the nucleon 
polarization terms in the response function (see Appendix E), and hence they can be evaluated to lowest order in the 
inelasticity. In the sum over iVp, In c .n p In c -i.n p will then give zero. For relativistic electrons we then find 



L^N^(N e > 1) = ~ 



iai; 



l N e ,N, 



,(wi/2eB) 1 



Pe,z \ 



Cy + 3c 2 A + 2c a (ca + cy)s n fl ■ B — 2c a (ca — cy)s p n ■ B 



Pe,z 

|A| 



T 2 



A-i,JV>x/2eS) 
Summing this expression over all N p gives 

L fiV N^(N e >l) = c v 

2c A (ca — Cy)s p fl ■ B 



{w\/2eB) 1 



Pe,z 

|A| 



(c 2 v - c 2 A )n ■ B. 



(D17) 



E 

N„=0 



3c A + 2ca(ca + cy)s n £l ■ B 



|A| 



(4-4)n-B. 



(D18) 
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In addition, it will be necessary to sum the matrix element against nucleon recoil terms from the response function 
which contains 



The needed summation rule is 



q 2 ± = eB(2N p + l)-p 2 ri 



J2SI 2 NS (x)=N- 



In our case this gives 



£ N P I 2 NeNp (w 2 ± /2eB) = N e + w 2 j2eB, 



(D19) 



(D20) 



(D21) 



JV„=o 



so that 



E 9x4.JV,K/2eB) = 2eB (N e + w 2 j2eB) + eB-p 2 n>x 
= eB(2N e + 1) + k\ + 2k ± ■ p n ± . 



N„=0 



(D22) 



Note that the large neutron momentum terms (p 2 ± ) cancelled so that the "averaged" recoil momentum has the 
expected size. The final result needed is then 



N p =0 

x 



2eBN e + k\+ 2k ± • p n , ± + eB^ 



E q±L^N^{N e > 1) 

\cy + 3c 2 A + 2c a (ca + c v )s n Vt ■ B - 2c a (ca - c v )s p fl 



B 



eB + (2eBN e + k\ + 2k ± • p n , ± ) 



(c 2 v - 4)n ■ b. 



(D23) 



APPENDIX E: NUCLEON RESPONSE FUNCTION FOR ABSORPTION 



The response function for absorption is defined in eq. (5J5). Using p PiZ to integrate over the z-momentum delta 



function gives p PtZ — p n . z + k z — p e ,z- Using p n ^ z to integrate over the energy delta function gives 

Qa + Q - q 2 /2m + (n Bp s p - HB n s n )B 

Pn.z = : , (>lj 

q z /m 

where we have defined the energy transfer by go = k — E e and the momentum transfer by q z = k z — p e , z , q\ = 
eB(2N p + 1) - Pnx, and q 2 = q\ + q 2 z . The result is 



/?? 



S SnSp = r-Jn(E n ) [1 - f p (E p )] , (E2) 



where the neutron energy is 



„ _Pl.± , Pn,z D Pl,± , [qo+Q-q 2 /2m+(nBpS p - HBnS n )B] 2 

n ~ ~2m~ 2m~~ ^ BnBSn = ^m~ + MqU^m) ^BnSnB, (E3) 

and the proton energy is E p = Q + E n + q^. Defining dimensionless parameters y, z via 



_ E n -H n _ fl n P„,_L [qo + Q~q 2 /2m + (nBpSp- fl B nS n )B] /J,BnS n B 

V ~ T ~ T 2mT 4T{q1/2m) T ' [ ' 
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and 



^ +9 ° + Q (E5) 



the response function can be written as 



The above expression is exact. We now consider the regime where the nucleons are nondegenerate. This is valid for 
the outer layers of the neutron star where asymmetric flux can develop. Since the nucleon spin energies are small, we 
shall expand S SnSp to first order in fisB. Using the nondegenerate nucleon conditions / P <1 and e~ y <C 1, we find: 

777 

S SnSp = ^e~y° (1 - Sy) , (E7) 



where y ~ yo + Sy, and 



yo 



Hn , Pl,± , (qo + Q- q 2 /2mf 



T 2mT AT(q 2 z /2m) 



HBnS n B ( q + Q-q 2 _/2m\ HB P s p B{ q + Q-q±/2m\ 
Sy = -^^{ 1+ ^ )'^~V qjj2^ J' (E8) 

The Sy term now contains all the dependence on the nucleon spins. This term will give rise to the nucleon contribution 
to the asymmetric parity violation effect. 



In evaluating the absorption opacity (eq. 5.8), it will be necessary to expand S SnS for small "inelasticity". If the 
nucleon mass were infinite, energy conservation would give exactly E e = k + Q. Thus we expect S SnS to be sharply 
peaked about this electron energy, with a width proportional to (T '/to) 1 / 2 . We can expand the electron energy around 
the peak in a series in the small parameter (T/m) 1 / 2 . There are two cases to consider: for the electron ground state 
(N e = 0), we shall want to define the dimensionless electron energy in terms of p e , z ] but for the case in which we are 
summing over a continuum of electron Landau levels it will be more convenient to define the dimensionless electron 
energy in terms of the perpendicular momentum p 2 . j_ = 2eBN e . 

In the N e — case, E e ~ \p e , z \ (neglecting to e ), we define the dimensionless electron energy u by 

p e , 2 =±(* + Q)(l + ctt), (E9) 



where 



2TV /2 Koj 



ml k + Q 



q z ,o = k z T(k + Q). (E10) 



(since e -C 1, we can set |1 + eu\ — 1 + eu over the interesting range of u). As discussed in §V.B, we will only need the 
dominant term for the N e = response function. Thus we can drop the nucleon polarization terms in Sy and work to 
the lowest order in inelasticity. With these approximations, we find 

^ P = r^ ex P[ ^-|b- u2 ) (^ = state), (Ell) 

which takes the form of a simple Gaussian in u. Since e -C 1, we can consider the range of u to extend from — oo to 
oo in the phase space integrals wi th e xponentially small error. 




When summing over N e in eq. (5.8), it will be more convenient to define u in terms of p\ j_ = 2eBN e . Let 

pl ± =2eBN e = El{l + eu), (E12) 

where 

Ei = (k+Qf - A ,„ e =(^)" J M!±« (E13 , 
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The allowed range of p e ^ z is now p e ^ z G [—(k + Q),(k + Q)] in order to keep E\ a positive number. The electron 
energy is 



The expressions needed for eq. are then 



1 



euE\ 



e 2 u 2 Ej_ 



eM-Vo) - exp - - — - u 



2(k + Q) 2 8(fc + Q) 4 

eu 3 El 2u(k + Q)q 2 ~* 
h 2(k + Q) 2 ~ 



emEl 



and 



8y = - 



x _q\_ eumEj | 



<& (k + Q)q 2 V 4(fc + Q) 2 

i + 4 + ^fi ^ 



2T 

2T ^ g 2 ' (fc + ^ 4 (fc + Q) 2 / 

Note that, as in the scattering case, there is a term in Sy with a coefficient which scales as 

emEl (mT) 1 / 2 /m\ 1 / 2 
(k + Q)q 2 ~ 

for k ~ T, and this term can be much greater than unity. 



/my 

It/ 



(E14) 



(E15) 



(E16) 



(E17) 



APPENDIX F: REPLACING THE SUMS OVER N E WITH INTEGRALS 



In this Appendix, we turn the sum over N e in eq. (5.8) into an integral. Let the sum be called 

oo 
JV e =l 

where the function T 1 is given by 

00 j2 



w)=d-/ e )E rmt e 

N P =0 J y ' s„,a p =±l 

Integrating the identity (see Ref. ]8Q|| ) 

OO OO 

(5(x — n) = exp(27rifca;) 

n= — oo fc= — oo 

against -F(a;) over the region x G [l,oo) gives the result 

00 ^ />oo 

^2 F(N e ) = --F(l) + dN e F(N e ) + oscillatory terms. 



(Fl) 



(F2) 



(F3) 



(F4) 



As discussed at the beginning of §|y|, we shall ignore the oscillatory terms. 

The F(l) terms will give expressions smaller than the integral over N e for both eq. ( 5.23| ) and (5.25). The reason 
is that the integral over N e effectively divides by eB, so that the factor of eB is eq. ( |5~8 ) is cancelled. For eq. (5.24), 
the F(l) term is odd in p e ^ z so that it integrates to zero (to lowest order in inelasticity) in the p e , z integral. The 
corrections involving inelasticity will make this term smaller by a factor of T/m than the N e — term. Our result is 
then 



1 = 



dN e F(N e ). 



(F5) 



When changing the variable of integration from N e to u using the expression 2eBN e — Ej_(l + eu) in cq.( 5.26| ), we 
can let u range from —00 to 00 since e ~ (T/m) 1 / 2 <C 1. 
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